Rec'dF6T/rT#3MAR 2005 



(12) INTERNATIONAL APPLICATION PUBLISHED UNDER THE PATENT COOPERATION TREA* 

(19) World Intellectual Property 
Organization 
international Bureau 

(43) International Publication Date 
2 June 2005 (02.06.2005) 




J TREATY (PCT) 

1CK528894 

lllIillllllllOIIIIllIlllDllllllIllllli 



PCT 



(10) International Publication Number 

WO 2005/049663 A2 



(51) International Pntcnt Classification 7 : C08F 10/00 

2/00, 2/34 

(21) International Application Number: 

PCT/US2003/014565 



8 May 2003 (08.05.2003) 
English 
English 



(22) International Filing Date: 

(25) Filing Language: 

(26) Publication Language: 

(30) Priority Data: 

PCT/US02/32767 15 October 2002 (15.10.2002) US 

(63) Related by continuation (CON) or continuation-in-part 
(CIP) to earlier application: 

U S PCT/US02/32767 (CIP) 

Piled on 15 October 2002 (15.10.2002) 

(71) Applicant (for all designated States except US): EXXON- 
MOBIL CHEMICAL PATENTS INC. [USAJS1; 5200 
Bay way Drive, Bay town, TX 77520-5200 (US). 



(72) Inventors; and 

(75) Inventors/Applicants (for US only): LONG, Robert, L. 
[US/US]; 15014 Cobre Valley, Houston, TX 77062 (US). 
1MPELMAN, Ryan, W. [US/USJ; 5200 Bayway Drive, 
Baytown, TX 77520 (US). CHANG, Shih, Y. [— /US]; 
15719 El Dorado Oaks, Houston, TX 77059 (US). AN- 
DREWS, Timothy, J. [US/US]; 5200 Bayway Drive, Bay- 
town, TX 77520 (US). YAHN, David, A. [US/US]; 20015 
Magnolia Rend, Humble, TX 77346 (US). MORROW, 
David [US/US]; 1015 Willowvale Drive, Taylor Lake Vil- 
lage, TX 77586 (US). 

(74) Agent: GRIFFIS, Andrew, B.; ExxonMobil Chemical 
Company, Law Technology, P.O. Box 2149, Baytown, TX 
77522-2149 (US). 

(81) Designated States (national): AE, AG, AL, AM, AT, AU, 
AZ, BA. BB, BG, BR, BY, BZ, CA, CH, CN, CO, CR, CU, 
CZ, DE, DK, DM, DZ, EC, EE, ES, FI, GB. GD, GE, GH, 
GM, HR, HU, ID, IL, IN, IS, JP, KB, KG, KP, KR, KZ, LC, 
LK, LR, LS, LT, LU, LV, MA, MD, MG, MK, MN, MW, 

[ Continued on next page] 



= (54 > Title; ON-LINE MEASUREMENT AND CONTROL OF POLYMER PROPERTIES BY RAMAN SPECTROSCOPY 

(57) Abstract: Methods are provided for 
determining and controlling polymer properties 
on-line in a polymerization reactor system, such 
as a fluidized bed reactor. The methods include 
obtaining a regression model for determining a 
polymer property, the regression model including 
principal component loadings and principal 
component scores, acquiring a Raman spectrum 
of a polyolefin sample comprising polyolefin, 
calculating a new principal component score from 
at least a portion of the Raman spectrum and the 
principal component loadings, and calculating the 
polymer pmnerLy by applying the new principal 
component score to the regression model. The 
property can be controlled by adjusting at least 
one polymerization parameter based on the 
calculated polymer property. 




0.916 D.918 0.920 

Lab D easily (g/an 3 ) 

JUrom density B 357 A 350 Linear (Ramon donaty^] 



0.922 



rt 

o 



O 

is 



0.925 
~ 0.923 -i 



| 0.921 
r 

E 0.919 

2 

0.S17 



0.915 




0.917 0.919 0.92J 0.923 . 

Lab demit} (g/cm 3 ) 



O Raman danary P 1107 A 6 101 X 6100 X 1003 Linear (Raman da»iry)| 




WO 2005/049663 A2 



MX, MZ, NO, NZ, OM, PH, PL, PT, RO, RU, SC, SD SB 
SG, SK, SU TJ, TM, TN, TR, TT. TZ, UA. UG. US, UZ,' 
VC, VN, YU. ZA. ZM, ZW. 

(84) Designated States (regional): ARIPO patent (GH, GM, 
KK, US, MW, MZ ; SI), SU SZ, TZ, UG, ZM, ZW), 
Eurasian patent (AM, AZ, BY, KG, KZ, MD, RU, TJ, TM), 
European patent (AT, BE, BG, CH, CY, CZ, DE, DK, EE, 
ES, FI, FR, GB, GR, HU, IE, IT, LU, MC. NL, PT, RO,' 




llllIBlllllllillliniBlllllMIll 

SE, SI, SK, TR), OAPI patent (BF, BJ, CF, CG, CI. CM, 
GA, GN, GQ, GW, ML, MR, NE. SN, TD, TG). 

Published: 

— without international search report and to be republished 
upon receipt of that report 

For two-letter codes and ot fie r abbreviations, refer to tlie "Guid- 
ance Notes on Codes and Abbreviations" appearing at the begin- 
ning of each regular issue of the PCI ' Gazette, 



WO 2005/049663 



1 



PCT/US2003/014565 



ON-LINE MEASUREMENT A ND CONTROL OF POLYMER PROPERTY 
BY RAMAN SPECTROSCOPY 

[0001] This application is a continuation-in-part and claims the benefit of 
5 international patent application PCT/US02/32767, filed October 15, 2002, which 
claims the benefit of U.S. Provisional Application No. 60/345,337, filed 
November 9, 2001. 

FIELD OF THE INVENTION 

10 [0002] The present invention is directed generally to methods of measuring 
polymer properties on-line in a polymerization reactor system, and using those 
measured properties to control the polymerization reaction. In particular, the 
present invention provides methods of measuring properties of polyolefins such as 
melt index and density on-line, using Raman spectroscopy, and methods of 

15 controlling a reactor using real-time, on-line polymer property data provided by 
Raman spectroscopic measurements. 

BACKGROUND 

[0003] Gas phase processes for the homopolymerization and copolymerization of 
20 monomers, especially olefin monomers, are well known in the art. Such processes 
can be conducted, for example, by introducing the gaseous monomer or 
monomers into a stirred and/or fluidized bed of resin particles and catalyst. 
[0004] In the fluidized-bed polymerization of olefins, the polymerization is 
conducted in a fluidized-bed reactor, wherein a bed of polymer particles is 
25 maintained in a fluidized state by means of an ascending gas stream including 
gaseous reaction monomer. The polymerization of olefins in a stirred-bed reactor 
differs from polymerization in a gas fluidized-bed reactor by the action of a 
mechanical stirrer within the reaction zone, which contributes to fluidization of 
the bed. As used herein, the term "fluidized-bed" also includes stirred-bed 
30 processes and reactors. 

[0005] The start-up of a fluidized bed reactor generally uses a bed of pre-formed 
polymer particles. During the course of polymerization, fresh polymer is 
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generated by the catalytic polymerization of the monomer, and polymer product is 
withdrawn to maintain the bed at constant volume. An industrially favored 
process employs a fluidization grid to distribute the fluidizing gas to the bed, and 
also to act as a support for the bed when the supply of gas is cut off. The polymer 
5 produced is generally withdrawn from the reactor via one or more discharge 
conduits disposed in the lower portion of the reactor, near the fluidization grid. 
The fluidized bed includes a bed of growing polymer particles, polymer product 
particles and catalyst particles. This reaction mixture is maintained in a fluidized 
condition by the continuous upward flow from the base of the reactor of a 
10 fluidizing gas which includes recycle gas drawn from the top of the reactor, 
together with added make-up monomer. 

[0006] The fluidizing gas enters the bottom of the reactor and is passed, 
preferably through a fluidization grid, upwardly through the fluidized bed. 
[0007] The polymerization of olefins is an exothermic reaction, and it is therefore 

15 necessary to cool the bed to remove the heat of polymerization. In the absence of 
such cooling, the bed would increase in temperature until, for example, the 
catalyst became inactive or the polymer particles melted and began to fuse. 
[0008] In the fluidized-bed polymerization of olefins, a typical method for 
removing the heat of polymerization is by passing a cooling gas, such as the 

20 fluidizing gas, which is at a temperature lower than the desired polymerization 
temperature, through the fluidized-bed to conduct away the heat of 
polymerization. The gas is removed from the reactor, cooled by passage through 
an external heat exchanger and then recycled to the bed. 

[0009] The temperature of the recycle gas can be adjusted in the heat exchanger to 
25 maintain the fluidized-bed at the desired polymerization temperature. In this 
method of polymerizing alpha olefins, the recycle gas generally includes one or 
more monomeric olefins, optionally together with, for example, an inert diluent 
gas or a gaseous chain transfer agent such as hydrogen. The recycle gas thus 
serves to supply monomer to the bed to fluidize the bed and to maintain the bed 
30 within a desired temperature range. Monomers consumed by conversion into 
polymer in the course of the polymerization reaction are normally replaced by 
adding make-up monomer to the recycle gas stream. 
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[0010] The material exiting the reactor includes the polyolefm and a recycle 
stream containing unreacted monomer gases. Following polymerization, the 
polymer is recovered. If desired, the recycle stream can be compressed and 
cooled, and mixed with feed components, whereupon a gas phase and a liquid 
phase are then returned to the reactor. 

[001 1J The polymerization process can use Ziegler-Natta and/or metallocene 
catalysts. A variety of gas phase polymerization processes are known. For 
example, the recycle stream can be cooled to a temperature below the dew point, 
resulting in condensing a portion of the recycle stream, as described in U.S. Patent 
Nos. 4,543,399 and 4,588,790. This intentional introduction of a liquid into a 
recycle stream or reactor during the process is referred to generally as a 
"condensed mode" operation. 

[0012] Further details of fluidized bed reactors and their operation are disclosed 
in, for example, U.S. Patent Nos. 4.243,619, 4,543,399, 5,352,749, 5,436,304, 
5,405,922, 5,462,999, and 6,218,484, the disclosures of which are incorporated 
herein by reference. 

[0013] The properties of the polymer produced in the reactor are affected by a 
variety of operating parameters, such as temperatures, monomer feed rates, 
catalyst feed rates, and hydrogen gas concentration. In order to produce polymer 

20 having a desired set of properties, such as melt index and density, polymer exiting 
the reactor is sampled and laboratory measurements carried out to characterize the 
polymer. If it is discovered that one or more polymer properties are outside a 
desired range, polymerization conditions can be adjusted, and the polymer 
resampled. This periodic sampling, testing and adjusting, however, is undesirably 

25 slow, since sampling and laboratory testing of polymer properties such as melt 
index, molecular weight distribution and density is time-consuming. As a result, 
conventional processes can produce large quantities of "off-spec" polymer before 
manual testing and control can effectively adjust the polymerization conditions. 
This occurs during production of a particular grade of resin as well as during the 

30 transition process between grades. 

[0014] Methods have been developed to attempt to provide rapid assessment of 
certain polymer properties and rapid adjustment of polymerization conditions. 
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PCT publications WO 01/09201 and WO 01/09203 disclose Raman-based 
methods using principal components analysis (PCA) and partial least squares 
(PLS) to determine concentrations of components in a slurry reactor. The 
concentration of a particular component, such as ethylene or hexene, is 
5 determined based on measurements of a known Raman peak corresponding to the 
component. U.S. Patent No. 5,999,255 discloses a method for measuring a 
physical property of a polymer sample, preferably nylon, by measuring a portion 
of a Raman spectrum of the polymer sample, determining a value of a preselected 
spectral feature from the Raman spectrum, and comparing the determined value to 
10 reference values. This method relies on identification and monitoring of 
preselected spectral features corresponding to identified functional groups, such as 
NH or methyl, of the polymer. 

[001 5] Additional background information can be found in U.S. Patent Nos. 
6,144,897 and 5,151,474; European Patent application EP 0 561 078; PCT 

15 publication WO 98/08066; and Ardell, G.G. et al, "Model Prediction for Reactor 
Control," Chemical Engineering Progress, American Institute of Chemical 
Engineers, U.S., vol. 79, no. 6, June 1, 1983, pages 77-83 (ISSN 0360-7275). 
[0016] It would be desirable to have methods of determining polymer properties 
such as melt index, density and molecular weight distribution, on-line in a 

20 fluidized bed polymerization reactor, without the need to preselect or identify 
spectral features of a polymer to monitor. It would also be desirable to have 
methods of controlling a gas-phase fluidized bed reactor to maintain desired 
polymer properties, based on a rapid, on-line determination of the polymer 
properties. 

25 

SUMMARY OF THE INVENTION 

[0017] In one aspect, the present invention provides a process for determining 
polymer properties in a polymerization reactor system. The process includes 
obtaining a regression model for determining a polymer property, the regression 
30 model including principal component loadings and principal component scores, 
acquiring a Raman spectrum of a polyolefm sample comprising polyolefin, 
calculating a new principal component score from at least a portion of the Raman 
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spectrum and the principal component loadings, and calculating the polymer 
property by applying the new principal component score to the regression model. 
10018] In another aspect, the present invention provides a process for controlling 
polymer properties in a polymerization reactor system. The process includes 
5 obtaining a regression model for determining a polymer property, the regression 
model including principal component loadings and principal component scores, 
acquiring a Raman spectrum of a polyolefin sample comprising polyolefin, 
calculating a new principal component score from at least a portion of the Raman 
spectrum and the principal component loadings, calculating the polymer property 
10 by applying the new principal component score to the regression model, and 
adjusting at least one polymerization parameter based on the calculated polymer 
property. In particular embodiments, the at least one polymerization parameter can 
be, for example, monomer feed rate, comonomer feed rate, catalyst feed rate, 
hydrogen gas feed rate, or reaction temperature. 
15 [0019] In one embodiment, the regression model is constructed by obtaining a 
plurality of Raman spectra of polyolefin samples, calculating principal component 
loadings and principal component scores from the spectra using principal 
component analysis (PCA), and forming the regression model using the principal 
component scores such that the regression model correlates the polymer property 
20 to the principal component scores. 

[0020] In another embodiment, the regression model is a locally weighted 
regression model. 

[0021] In another embodiment, the method includes: obtaining a first regression 
model for determining a first polymer property, the first regression model 

25 including first principal component loadings and first principal component scores; 
obtaining a second regression model for determining a second polymer property, 
the second regression model including second principal component loadings and 
second principal component scores; acquiring a Raman spectrum of a sample 
comprising polyolefin; calculating a new first principal component score from at 

30 least a portion of the Raman spectrum and the first principal component loadings; 
calculating a new second principal component score from at least a portion of the 
Raman spectrum and the second principal component loadings; calculating the 
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first polymer property by applying the new first principal component score to the 
first regression model; and calculating the second polymer property by applying 
the new second principal component score to the second regression model. 
[0022] In another embodiment, the sample includes polyolefin particles. 
5 [0023] In another embodiment, the Raman spectrum is acquired by providing a 
sample of polyolefin particles and irradiating the sample and collecting scattered 
radiation during a sampling interval using a sampling probe, wherein there is 
relative motion between the sample and the sampling probe during at least a 
portion of the sampling interval. The relative motion serves to effectively 

10 increase the field of view of the sampling probe, providing more accurate data. 

[0024] In another embodiment, the Raman spectrum is acquired from a probe 
inserted into the reactor or downstream from the reactor. In a preferred 
embodiment the reactor is a gas phase polymerization reactor and more preferably 
is a fluidized bed reactor, e.g., a Unipol reactor or a gas phase, fluidized bed 

1 5 reactor having an optional cyclone. 

[0025) In other embodiments, suitable polymer properties include, for example, 
density, melt flow rates such as melt index or flow index, molecular weight, 
molecular weight distribution, and various functions of such properties. 

20 BRIEF DESCRIPTION OF THE DRAWINGS 

[0026] Figure 1 is a block diagram of a gas-phase reactor. 

[0027] Figure 2 is a block diagram of a Raman analyzer according to the 
invention. 

[0028] Figure 3 illustrates one embodiment of a fiber optic Raman probe. 
25 [0029] Figure 4 illustrates one embodiment of a sample chamber. 

[0030] Figure 5 is a representative Raman spectrum of a granular linear low 
density polyethylene polymer sample. 

[0031] Figures 6a and 6b show predicted versus measured melt indices in low and 
high melt index ranges, respectively, according to Examples 1 and 2. 
30 [0032] Figure 7 shows predicted versus measured density according to 
Example 3. 
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10033] Figures 8a and 8b show predicted versus measured melt indices from on- 
line Raman analyses in metallocene- and Ziegler-Natta-catalyzed reactions, 
respectively, according to Examples 4-5. 

[0034] Figures 9a and 9b show predicted versus measured density from on-line 
Raman analyses in metallocene- and Ziegler-Natta-catalyzed reactions, 
respectively, according to Examples 6-7. 

[0035] Figure 10 shows predicted versus measured melt indices from on-line 
Raman analyses in a commercial-scale fluidized-bed reactor, over a period of 
about five weeks. 

[0036] Figure 1 1 shows predicted versus measured densities from on-line Raman 
analyses in a commercial-scale fluidized-bed reactor, over a period of about five 
weeks. 



DETAILED DF.S ("!R TPTTO'M 

[0037JIn one embodiment, the present invention provides a method of 
determining polyolefm polymer properties on-line, i.e., as the polyolefin is 
produced in a reactor system, without the need for external sampling and analysis. 
The method includes obtaining a regression model for determining a polymer 
property, the regression model including principal component loadings and 
principal component scores, acquiring a Raman spectrum of a polyolefin sample, 
calculating a new principal component score from at least a portion of the Raman' 
spectrum and the principal component loadings, and calculating the polymer 
property by applying the new principal component score to the regression model. 
[0038] In one embodiment, the method is used to determine polymer properties 
25 on-line in a fluidized-bed reactor system. Fluidized-bed reactors are well-known 
in the art; a particular, non-limiting example of a fluidized bed reactor is described 
herein, for illustrative purposes only. Those skilled in the art will recognize that 
numerous modifications and enhancements can be made, as desired, to the 
fluidized-bed reactor. 
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Fluidized-Bed Reactor 

[0039] FIG. 1 illustrates one example of a gas phase polymerization reactor 
system which comprises a gas-phase' fluidized bed reactor 20 having a reactor 
5 body 22, which is generally an upright cylinder having a fluidization grid 24 
located in its lower regions. The reactor body 22 encloses a fluidized bed zone 26 
and a velocity reduction zone 28 which is generally of increased diameter 
compared to the diameter of the fluidized bed zone 26 of the reactor body 22. 
[0040] The gaseous reaction mixture leaving the top of the reactor body 22, 

10 termed the "recycle gas stream," contains principally unreacted monomer, 
unreacted hydrogen gas, inert condensable gases such as isopentane, and inert 
non-condensable gases such as nitrogen. The recycle gas stream is transferred via 
line 30 to compressor 32, and from compressor 32 to heat exchanger 34. An 
optional cyclone separator 36 may be used as shown, preferably upstream of 

15 compressor 32, to remove fines, if desired. An optional gas analyzer 38 can be 
used if desired, to sample the recycle gas stream to determine concentrations of 
various components. Typically, the gas analyzer is a gas phase chromatograph 
(GPC), or a spectrograph such as a near-infrared spectrometer or a fourier 
transform near-infrared spectrometer (FT-NIR). An additional heat exchanger 

20 (not shown) may also be used if desired, preferably upstream of compressor 32. 
[0041] The cooled recycle gas stream exits the heat exchanger 34 via line 40. As 
discussed above, the cooled recycle gas stream can be gaseous, or can be a 
mixture of gaseous and liquid phases. Figure 1 shows an optional configuration 
wherein at least a portion of the recycle gas stream is cooled to a temperature at or 

25 below the temperature where liquid condensate begins to form (the dew point). 
All or a portion of the resultant gas liquid mixture is transferred via line 40 to a 
separator 42, where all or a portion of the liquid is removed. All or a portion of the 
gas stream, which may contain some liquid, is transferred via line 44 to a point 
below the fluidization grid 24 in the lower region of the reactor. An amount of 

30 upwardly flowing gas, sufficient to maintain the bed in a fluidized condition, is 
provided in this way. 
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[0042] Those skilled in the art will understand that less gas is required to maintain 
fluidization when the reactor employed is a stirred bed reactor. 
[0043] An optional compressor 46 may be provided to ensure that a sufficient 
velocity is imparted to the gases flowing through line 44 into the bottom of the 
reactor. The gas stream entering the bottom of the reactor may contain condensed 
liquid, if desired. 

[0044] All or a portion of the liquid phase separated from the recycle stream in 
separator 42 is transferred via line 48 to a manifold 50 located at or near the top of 
the reactor. If desired, a pump 52 may be provided in line 48 to facilitate the 
transfer of liquid to manifold 50. The liquid entering manifold 50 flows downward 
into manifold 54 through a plurality of conduits 56 which have good heat 
exchange properties and which are in heat exchange contact with the wall of the 
reactor. The passage of liquid through the conduits 56 cools the interior wall of 
the reactor and warms the liquid to a greater or lesser extent, depending upon the 
temperature differential and the duration and extent of heat exchange contact. 
Thus by the time the liquid entering manifold 50 reaches manifold 54, it has 
become a heated fluid which may have remained in an entirely liquid state or it 
may have become partially or totally vaporized. 

[0045] As shown in FIG. 1, the heated fluid (gas and/or liquid) is passed from 
manifold 54 via line 58 to combine with gases leaving the separator 42 via line 44, 
prior to entry into the reactor in the region below the fluidization grid 24. In like 
manner, make-up monomer can be introduced into the reactor in either liquid or 
gaseous form via line 60. Gas and/or liquid collected in manifold 54 may also be 
transferred directly into the reactor (not shown) in the region below the 
25 fluidization grid. 

[0046] Product polymer particles can be removed from the reactor via line 62 in 
the conventional way, as for example by the method and apparatus described in 
U.S. Pat. No. 4,621,952. Although only one line 62 is shown in the Figure, 
typical reactors can include more than one line 62. 
30 [0047] Catalyst is continuously or intermittently injected into the reactor using a 
catalyst feeder (not shown) such as the device disclosed in U.S. Pat. No. 
3,779,712. The catalyst is preferably fed into the reactor at a point 20 to 40 
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percent of the reactor diameter away from the reactor wall and at a height of about 
5 to about 30 percent of the height of the bed. The catalyst can be any catalyst 
suitable for use in a fluidized bed reactor and capable of polymerizing ethylene, 
such as one or more metallocene catalysts, one or more Ziegler-Natta catalysts,' 
bimetallic catalysts, or mixtures of catalysts. 

[0048] A gas which is inert to the catalyst, such as nitrogen or argon, is preferably 
used to carry catalyst into the bed. Cold condensed liquid from either separator 42 
or from manifold 54 may also be used to transport catalyst into the bed. 
10049] In methods of the present invention, the fluidized bed reactor is operated to 
form at least one polyolefin homopolymer or copolymer. Suitable polyolefins 
include, but are not limited to, polyethylene, polypropylene, polyisobutylene, and 
homopolymers and copolymers thereof. 

]0050] In one embodiment, the at least one polyolefin includes polyethylene 
homopolymer and/or copolymer. Low density polyethylene ("LDPE") can be 
prepared at high pressure using free radical initiators, or in gas phase processes 
using Ziegler-Natta or vanadium catalysts, and typically has a density in the range 
of 0.916-0.940 g/cm 3 . LDPE is also known as "branched" or "heterogeneously 
branched" polyethylene because of the relatively large number of long chain 
branches extending from the main polymer backbone. Polyethylene in the same 
density range, i.e., 0.916 to 0.940 g/cm 3 , which is linear and does not contain long 
chain branching is also known; this "linear low density polyethylene" ("LLDPE") 
can be produced with conventional Ziegler-Natta catalysts or with metallocene 
catalysts. Relatively higher density LDPE, typically in the range of 0.928 to 
0.940g/cm 3 , is sometimes referred to as medium density polyethylene ("MDPE"). 
25 Polyethylenes having still greater density are the high density polyethylenes 
("HDPEs"), i.e, polyethylenes having densities greater than 0.940 g/cm 3 , and are 
generally prepared with Ziegler-Natta catalysts. Very low density polyethylene 
("VLDPE") is also known. VLDPEs can be produced by a number of different 
processes yielding polymers with different properties, but can be generally 
30 , described as polyethylenes having a density less than 0.91 6 g/cm 3 , typically 0.890 
to 0.915 g/cm 3 or 0.900 to 0.915 g/cm 3 . 
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[0051] Polymers having more than two types of monomers, such as terpolymers, 
are also included within the term "copolymer" as used herein. Suitable 
comonomers include a-olefins, such as C 3 -C 20 a-olefins or C3-C12 a-olefins. The 
a-olefin comonomer can be linear or branched, and two or more comonomers can 
5 be used, if desired. Examples of suitable comonomers include linear C3-C12 
a-olefins, and a-olefins having one or more C1-C3 alkyl branches, or an aryl 
group. Specific examples include propylene; 3 -methyl- 1-butene; 
3,3-dimethyl-l-butene; 1-pentene; 1-pentene with one or more methyl, ethyl or 
propyl substituents; 1-hexene with one or more methyl, ethyl or propyl 

10 substituents; 1-heptene with one or more methyl, ethyl or propyl substituents; 
1-octene with one or more methyl, ethyl or propyl substituents; 1-nonene with one 
or more methyl, ethyl or propyl substituents; ethyl, methyl or dimethyl-substituted 
1-decene; 1-dodecene; and styrene. It should be appreciated that the list of 
comonomers above is merely exemplary, and is not intended to be limiting. 

15 Preferred comonomers include propylene, 1-butene, 1-pentene, 4-methyl-l- 
pentene, 1-hexene, 1-octene and styrene. 

[0052] Other useful comonomers include polar vinyl, conjugated and non- 
conjugated dienes, acetylene and aldehyde monomers, which can be included in 
minor amounts in terpolymer compositions. Non-conjugated dienes useful as co- 
20 monomers preferably are straight chain, hydrocarbon diolefins or cycloalkenyl- 
substituted alkenes, having 6 to 15 carbon atoms. Suitable non-conjugated dienes 
include, for example: (a) straight chain acyclic dienes, such as 1,4-hexadiene and 

1.6- octadiene; (b) branched chain acyclic dienes, such as 5-methyl- 1,4-hexadiene; 

3.7- dimethyl-l,6-octadiene; and 3,7-dimethyl-l,7-octadiene; (c) single ring 
25 alicyclic dienes, such as 1,4-cyclohexadiene; 1 ,5-cycIo-octadiene and 1,7- 

cyclododecadiene; (d) multi-ring alicyclic fused and bridged ring dienes, such as 
tetrahydroindene; norbornadiene; methyl-tetrahydroindene; dicyclopentadiene 
(DCPD); bicyclo-(2.2.1)-hepta-2,5-diene; alkenyl, alkylidene, cycloalkenyl and 
cycloalkylidene norbornenes, such as 5-methylene-2-norbornene (MNB), 5- 
30 propenyl-2-norbornene, 5-isopropylidene-2-norbomene, 5-(4-cyclopentenyl)-2- 
norbornene, 5-cyclohexylidene-2-norbornene, and 5-vinyl-2-norbornene (VNB); 
and (e) cycloalkenyl-substituted alkenes, such as vinyl cyclohexene, allyl 
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cyclohexene, vinyl cyclooctene, 4-vinyl cyclohexene, allyl cyclodecene, and vinyl 
cyclododecene. Of the non-conjugated dienes typically used, the preferred dienes 
are dicyclopentadiene, 1,4-hexadiene, 5-methylene-2-norbornene, 5-ethylidene-2- 
norbornene, and tetracyclo-(A-l l,12)-5,8-dodecene. Particularly preferred 
5 dioleflns are 5-ethyIidene-2-norbornene (ENB), 1,4-hexadiene, dicyclopentadiene 
(DCPD), norbornadiene, and 5-vinyl-2-norbornene (VNB). 
[0053] The amount of comonomer used will depend upon the desired density of 
the polyolefm and the specific comonomers selected. One skilled in the art can 
readily determine the appropriate comonomer content appropriate to produce a 
10 polyolefin having a desired density. 

Raman Spectroscopy 

[0054J Raman spectroscopy is a well-known analytical tool for molecular 
characterization, identification, and quantification. Raman spectroscopy makes 

1 5 use of inelastically scattered radiation from a non-resonant, non-ionizing radiation 
source, typically a visible or near-infrared radiation source such as a laser, to 
obtain information about molecular vibrational-rotational states. In general, non- 
ionizing, non-resonant radiation is scattered elastically and isotropically (Raleigh 
scattering) from a scattering center, such as a molecule. Subject to well-known 

20 symmetry and selection rules, a very small fraction of the incident radiation can be 
inelastically and isotropically scattered, with each inelastically scattered photon 
having an energy E = hn 0 ± | Ej-j. - Ey | , where hu 0 is the energy of the incident 
photon and | Brj - E y | is the absolute difference in energy between the final (i* j 1 ) 
and initial (ij) vibrational-rotational states of the molecule. This inelastically 

25 scattered radiation is the Raman scattering, and includes both Stokes scattering, 
where the scattered photon has lower energy than the incident photon (E = hn 0 - 
| Ej'j. - Ey | ), and anti-Stokes scattering, where the scattered photon has higher 
energy than the incident photon (E = hn 0 + | Eyj - E,j | ). 

[0055] Raman spectra are typically shown as plots of intensity (arbitrary units) 
30 versus "Raman shift", where the Raman shift is the difference in energy or 
wavelength between the excitation radiation and the scattered radiation. The 
Raman shift is typically reported in units of wavenumbers (cm" 1 ), i.e., the 



WO 2005/049663 ^ PCT/US2003/014565 



10 



15 



20 



25 



30 



reciprocal of the wavelength shift in centimeters. Energy difference 
lEpj.-Eijland wavenumbers (oo) are related by the expression | Ejj. - E, j | = 
hcco, where h is Planck's constant, c is the speed of light in cm/s, J „ is the 
reciprocal of the wavelength shift in centimeters.- - 

[0056] The spectral range of the Raman spectrum acquired is not particularly 
limited, but a useful range includes Raman shifts (Stokes and/or anti-Stokes) 
corresponding to a typical range of polyatomic vibrational frequencies, generally 
from about 100 cm" to about 4000 cm 1 . It should be appreciated that useful 
spectral information is present in lower and higher frequency regions. For 
example, numerous low frequency molecular modes contribute to Raman 
scattering in the region below 100 cm ' Raman shift, and overtone vibrations 
(harmonics) contribute to Raman scattering in the region above 4000 cm" 1 Raman 
shift. Thus, if desired, acquisition and use of a Raman spectrum as described 
herein can include these lower and higher frequency spectral regions. 
[0057] Conversely, the spectral region acquired can be less than all of the 100 
cm' 1 to 4000 cm 1 region. For many polyolefins, the majority of Raman scattering 
intensity will be present in a region from about 500 cm" 1 to about 3500 cm' 1 or 
from 1000 cm' 1 to 3000 cm". The region acquired can also include a plurality of 
sub-regions that need not be contiguous. 

[0058] As explained below, it is a particular advantage of the methods described 
herein that Raman scattering intensity data is useful in determining properties of 
polyolefin particles without the need to identify, select, or resolve particular 
spectral features. Thus, it is not necessary to identify a particular spectral feature 
as being due to a particular mode of a particular moiety of the polyolefin, nor is it 
necessary to selectively monitor Raman scattering corresponding to a selected 
spectral feature. Indeed, it has been surprisingly found that such selective 
monitoring disadvantageous^ disregards a wealth of information content 
embedded in the spectrum that, heretofore, has generally been considered to be 
merely unusable scattering intensity disposed between and underlying the 
identifiable (and thus presumed useful) bands. Accordingly, in the methods 
described herein, the Raman spectral data acquired and used includes a plurality 
of frequency or wavelength shift, scattering intensity (x, y) measurements over 
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relatively broad spectral regions, including regions conventionally identified as 
spectral bands and regions conventionally identified as interband, or unresolved 
regions. 

[0059J The frequency spacing of acquired data can be readily determined by one 
skilled in the art, based on considerations of machine resolution and capacity 
acquisition time, data analysis time, and information density. Similarly, the' 
amount of signal averaging used is readily determined by one skilled in the art 
based on machine and process efficiencies and limitations. 
[0060] The spectral region measured can include Stokes scattering (i.e., radiation 
scattered at frequencies lower than the excitation frequency), anti-Stokes 
scattering (i.e., radiation scattered at frequencies higher than the excitation 
frequency), or both. Optionally, polarization information embedded in the Raman 
scattering signal can also be used, and one skilled in the art readily understands 
how to acquire Raman polarization information. However, determining polymer 
properties as described herein does not require the use of polarization information. 
In some embodiments described herein, any Raman polarization is essentially 
randomized as a result of interactions with the fiber optic conduit used to convey 
the signal to the signal analyzer, as described below. 

20 Raman Instrumentation 

[0061] Referring now to Figure 2, the instrumentation used to collect and process 
Raman data includes a Raman subsystem 100, a sample subsystem 200, and a data 
subsystem 300. As shown in Figure 2, the sample subsystem 200 is in 
communication with reactor 20 via polymer output line 62 (see also Figure 1). 

25 Each of these subsystems is described below. Figure 2 illustrates the case wherein 
the sampling is extractive, as opposed to the case wherein the sampling is in situ, 
as described in more detail below. 

[0062] In an embodiment not shown, the Raman probe 204 may be inserted 
directly into reactor body 22. Reactor body 22 may thus act as sample subsystem 
30 200. It will be recognized by one of skill in the art in possession of the present 
disclosure that Raman probe 204 may be used anywhere in process where granular 
resin could be collected and analyzed by a Raman probe or anywhere in the 
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process where granular resin can move relative to a Raman probe, e.g., in the 
cycle gas piping (e.g., line 30 in Fig. 1), in the product discharge system 
downstream of the exiting point of product (e.g., from 22 into line(s) 62 in Fig. 1), 
in the transfer line between the product discharge system and the 
5 purger(s)/degasser(s), in one or more of the purger(s)/degasser(s), in the transfer 
line to finishing/pack-out, in the feed bins to the extruder/mixer (not shown), and 
the like. 

[0063] In an embodiment not shown, Raman probe 204 is inserted into fluidized 
bed zone 26, more preferably in the lower half of zone 26 but above grid 24. 

10 

Raman Subsystem 

[0064] The Raman subsystem includes a Raman spectrometer, the principal 
components of which are an excitation source 102, a monochromator 104, and a 
detector 106. Raman spectrometers are well-known analytical instruments, and 

1 5 thus only a brief description is provided herein. 

[0065] A Raman spectrometer includes an excitation source 102 which delivers 
excitation radiation to the sample subsystem 200. Scattered radiation is collected 
within the sample subsystem 200 (described below), filtered of Raleigh scattered 
light, and dispersed via monochromator 104. The dispersed Raman scattered light 

20 is then imaged onto a detector 106 and subsequently processed in data subsystem 
300, as further described below. 

Excitation Source 

[0066] The excitation source and frequency can be readily determined based on 
25 considerations well-known in the art. Typically, the excitation source 102 is a 
visible or near infrared laser, such as a frequency-doubled NdrYAG laser (532 
nm), a helium-neon laser (633 nm), or a solid-state diode laser (such as 785 nm). 
The laser can be pulsed or continuous wave (CW), polarized as desired or 
randomly polarized, and preferably single-mode. Typical excitation lasers will 
30 have 100 to 400 mW power (CW), although lower or higher power can be used as 
desired. Light sources other than lasers can be used, and wavelengths and laser 
types and parameters other than those listed above can also be used. It is well- 
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known that scattering, including Raman scattering, is proportional to the fourth 
power of the excitation frequency, subject to the practical limitation that 
fluorescence typically overwhelms the relatively weak Raman signal at higher 
frequencies. Thus, higher frequency (shorter wavelength) sources are preferred to 
5 maximize signal, while lower frequency (longer wavelength) sources are preferred 
to minimize fluorescence. One skilled in the art can readily determine the 
appropriate excitation source based on these and other considerations, such as 
mode stability, maintenance time and costs, capital costs, and other factors well 
understood in the art. 

10 [0067] The excitation radiation can be delivered to the sample subsystem 200, and 
the scattered radiation collected from the sample subsystem, by any convenient 
means known in the art, such as conventional beam manipulation optics, or fiber 
optic cables. For an on-line process measurement, it is particularly convenient to 
deliver the excitation radiation and collect the scattered radiation fiber-optically. 

15 It is a particular advantage of Raman spectroscopy that the excitation radiation 
typically used is readily manipulated fiber optically, and thus the excitation source 
can be positioned remotely from the sampling region. A particular fiber optic 
probe is described below; however, one skilled in the art will appreciate that the 
Raman system is not limited to any particular means of radiation manipulation. 

20 

Monochromator 

[0068] The scattered radiation is collected and dispersed by any convenient means 
known in the art, such as a fiber optic probe as described below. The collected 
scattered radiation is filtered to remove Raleigh scattering and optionally filtered 

25 to remove fluorescence, then frequency (wavelength) dispersed using a suitable 
dispersive element, such as a blazed grating or a holographic grating, or 
interferometrically (e.g., using Fourier transforms). The grating can be fixed or 
scanning, depending upon the type of detector used. The monochromator 104 can 
be any such dispersive element, along with associated filters and beam 

30 manipulation optics. 
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[0069] The dispersed Raman scattering is imaged onto a detector 106. The choice 
of detector is easily made by one skilled in the art, taking into account various 
factors such as resolution, sensitivity to the appropriate frequency range, response 
time, etc. Typical detectors include array detectors generally used with fixed- 
dispersive monochromators, such as diode arrays or charge coupled devices 
(CCDs), or single element detectors generally used with scanning-dispersive 
monochromators, such as lead sulfide detectors and indium-gallium-arsenide 
detectors. In the case of array detectors, the detector is calibrated such that the 
frequency (wavelength) corresponding to each detector element is known. The 
detector response is delivered to the data subsystem 300 which generates a set of 
frequency shift, intensity (x,y) data points which constitute the Raman spectrum. 

Sample Subsystem 

[0070] The sample subsystem 200 couples the Raman subsystem 100 to the 
polymerization process. Thus, the sample subsystem 200 delivers the excitation 
radiation from the excitation source 102 to the polymer sample, collects the 
scattered radiation, and delivers the scattered radiation to the monochromator 104. 
[0071] As noted above, the excitation radiation can be delivered to and collected 
from the polymer sample by any convenient means, such as using conventional 
optics or fiber optic cables. 

[0072] In one embodiment, the sample subsystem includes a probe 204 and a 
sample chamber 202. Figure 3 shows a block diagram of one embodiment of a 
fiber optic probe. The probe includes a fiber optic bundle 206 including one or 
more fiber optic cables 208 carrying the excitation radiation from the excitation 
source toward the sample, and one or more fiber optic cables 210 carrying the 
collected scattered radiation from the sample. Fiber optic cables 208 are in optical 
communication with the excitation source (102 in Figure 2), and fiber optic cables 
210 are in optical communication with the monochromator (104 in Figure 2). The 
excitation and scattered radiation can be manipulated using well-known 
techniques. Thus, it should be appreciated that the particular optical setup shown 
in Figure 3 is merely exemplary. Excitation radiation 212 is directed via optics 
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214 to a holographic grating 216 and spatial filter 218 to remove silica Raman due 
to the fiber optic cable, then directed via mirror 220 and beam combiner 222 to 
sampling optics 224 and sample chamber 202. Scattered radiation is collected via 
sampling optics 224 and directed through beam combiner 222, a notch filter 226 
5 to remove the Raleigh scattered radiation, and into fiber optic cables 210. 

[0073] The sample in the sample chamber includes a plurality of polymer 
particles (granules), and represents the polymer product as discharged from the 
reactor. Advantageously, it is not necessary that the sample be free of liquid- 
phase components, such as residual solvent or other liquid hydrocarbons that may 

1 0 be present in the polymer in the discharge line of a fluidized-bed reactor. 

[0074] Raman probes such as described herein are imaging, in that they have a 
focused field of view. An imaging probe is the most efficient optical 
configuration, and because the Raman signal is weak the imaging probe collects 
as much scattered light as possible. A disadvantage of an imaging probe is that 

15 the probe "sees" only a very small amount of the sample at any one time. For a 
typical fluidized-bed process, a fixed imaging probe has a field of view 
corresponding to only 1 or 2 polymer granules. Thus, the data collected in a static 
mode may not be representative of the bulk material. 

[0075] In one embodiment, the disadvantage of a limited field of view is 
20 overcome by providing relative motion between the sample and the Raman probe, 
so that the probe collects scattering from many polymer granules over the course 
of the sampling interval. Thus, for example, the probe can be moved through the 
sample during at least a portion of the sampling interval or, equivalently, the 
sample or sample chamber can be moved relative to a fixed probe during at least a 
25 portion of the sampling interval, or both can be moved. In a particular 
embodiment, it is convenient to keep the sample chamber stationary and move the 
Raman probe into and out of the sample chamber during the sampling interval by 
linearly translating the probe using a linear actuator. One skilled in the art will 
readily appreciate, however, that relative motion between the sample granules and 
30 the probe can be achieved by numerous other mechanisms, such as, for example, 
allowing polymer granules to pass by a stationary probe, as would occur, for 
example, if the Raman probe is inserted within the reactor body 22. Additional 
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embodiments can thus be envisioned wherein the Raman probe is placed or 
inserted in situ into the polymerization reactor system where granular polymer is 
moving, i.e., without the need for a sample chamber, such as shown by sample 
chamber 202 in Fig. 2, being added to the system. As used herein the term 
"polymerization reactor system" encompasses the system illustrated (but not 
limited) by Fig. 1, that is, from the polymerization reactor body (20 in Fig. 1) 
through the finishing apparatus prior to extruding/pelletization. It is preferred that 
the polymerization reactor system be a gas phase polymerization reactor system. 
[0076] In an embodiment, particularly in the case where polymer granules pass by 
a probe, whether or not the probe is held stationary, the probe may be purged of 
collected polymer product to prevent the field of view from getting coated. This 
may be accomplished, for instance, by the use of a purge with N 2) H 2 , ethylene, 
isopentane, hexane, mineral oil, n-butane and the like. In an even more preferred 
embodiment, there is a cycling between periods of data collection and probe purge 
in order to obtain optimal readings. The probe purge/data collection cycle times 
may be varied, and the depth of the probe insertion may also be varied. 
[0077] One advantage of inserting the probe directly into reactor body 22 is 
earlier indication of polymer properties and also problems with reactor 
operability, such as onset of sheeting or fouling which may cause the probe tip to 
20 plug. 

[0078] Appropriate probes are available commercially, for instance, from Axiom 
Analytical, Inc. and Kaiser Optical Systems, Inc. 

[0079] As a specific example, a particular sampling system used in Examples 4-7 
below is now described. It should be appreciated that this particular system is 
25 exemplary and not limiting. 

[0080] A fluidized-bed polymerization plant having two reactors was used, with 
one reactor producing metallocene-catalyzed LLDPE resin, and the other reactor 
producing Ziegler-Natta catalyzed LLDPE resin. Referring now to Figure 4, each 
reactor 20 (only one reactor shown) has two dump valves A and B that alternate to 
remove product from the reactor. The product is pneumatically conveyed through 
product discharge pipe 62 with 90 psi (0.6 MPa) nitrogen at a speed of about 60 
miles per hour (0.4 m/s). At this speed the slug of product dumped from a reactor 
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will only be present at any one point in the pipe for a few seconds. However, it is 
preferred to average the Raman signal for 60-120 seconds' to improve the signal- 
to-noise ratio. To accomplish this, a small amount of product (about 800 grams) 
is trapped and helci in a sample chamber 202 as the slug passes through the 
5 product discharge pipe 62. The sample chamber 202 is attached to the product 
discharge pipe 62 by a 1 inch (25 mm) diameter pipe 62b and a pneumatically 
actuated valve C or D. The operation of the valves C and D is controlled by the 
Raman analyzer, but could also be controlled by an auxiliary system. The Raman 
analyzer waits for a signal from the reactor telling it that the dump valve A or B 
10 has opened. The Raman analyzer then opens valve C or D connecting the sample 
chamber 202 to the product discharge pipe 62, and waits for a time predetermined 
to be sufficient to have allowed the slug of product to have passed by the sample 
capture point. The Raman analyzer next closes the sample capture valve C or D, 
trapping the captured sample of product in the sample chamber 202. 
15 [0081] The Raman analyzer probe 204 includes a probe head 230 enclosing the 
filtering and optical (not electronic) signal processing elements, and a sample 
interface 232, which is an 8" long by 0.5" diameter (20 cm x 1.3 cm) tube. Tube 
232 is inserted through the end of the sample chamber opposite to where the 
sample enters, so that it comes in contact with the sample. A pneumatic linear 
20 actuator 234 is attached to the probe 204 to slowly draw the probe out of the 
sample chamber and then reinsert it during a sample collection interval. This 
probe movement causes sample to flow across the front of the probe, providing a 
continually changing sample for measurement. 

[0082] The reactor 20 dumps on a 3-6 minute cycle (grade dependent), alternating 
25 between 2 lines 62 controlled by valves A and B. Sample is collected from only 

one of the lines. The sample system operates by waiting for a Sample Ready 

signal from the reactor telling the Raman analyzer that a sample is being dumped. 

The Sample Ready signal is in the form of a digital input to the Raman analyzer. 

When the analyzer receives the Sample Ready signal, there is a sequence of tasks 
30 it performs prior to setting up the valves for the Capture Sample operation, which 

are: 
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[0083] Check to determine if the Sample Ready is for the next stream. In the 
Raman control software, there is a stream sequence list that the operator sets to 
tell the analyzer which reactor(s) to sample and measure. Typically, this would be 
1,2,1,2, etc., for a two reactor system, but under some circumstances such as a 
5 grade transition on reactor 1, the operator might want to sample, for example, 
1,1,1,2,1,1,1,2, etc. Thus, the analyzer checks to make sure the dump indicator it 
receives is consistent with the current stream sequence. If no t, the analyzer 
ignores the signal. 

[0084] Check that the Reactor On-line digital input for this reactor is valid. The 
10 typical stream sequence 1,2,1,2 ... may be in effect, but the operator may decide to 
only monitor a single reactor, such as during a transition or upset. The reactor 
receives separate digital inputs for each reactor, which tell it whether or not to 
sample a particular reactor regardless of the active or current stream sequence list. 
[0085] Wait a set time interval between the Sample Ready signal and setting 
15 valves for Capture Sample. 

[0086] Set Valves for Capture Sample. 

[0087] The valve states are shown in the table below for a sequence sampling 
through the A valve of product discharge line 62, with state "C" being closed, and 
state "O" being open. 



Valve States For Sampling 



Valve 


A 


B 


c 


D 


E 


F 


Waiting for Sample 


C 


C 


c 


C 


c 


C 


Capture Sample 


O 


C 


o 


O 


c 


C 


Measure Spectrum 


C 


C 


c 


c 


c 


O 


Eject Sample 


C 


C 


o 


c 


0 


O 


Reset Probe 


c 


c 


0 


c 


0 


C 



[0088] Sample Capture is accomplished by opening the sample chamber valves C 
and D. In the configuration where product is discharged through the A valve of 
25 product discharge line 62, an open valve C permits the sample to enter sample 
chamber 202, and an open valve D serves as a vent. A portion of the discharged 
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polymer product in 90 psig nitrogen being transported at about 60 miles per hour 
packs into the sample chamber 202 attached to a bend in the product discharge 
line 62. Once the sample chamber 202 is full, the analyzer performs a series of 
operations to complete data collection and prepare for the next sample. These 
5 operations include: 

[0089J Wait a specified time interval /after the Capture Sample valve state is set. 
[0090] Set the Measure Spectrum valve state. 
[0091] Eject the sample 
[0092] Reset the Probe Position. 
1 0 [0093] Set the Waiting for Sample valve state 
[0094] Update the stream sequence information. 

[0095] The probe is attached to linear actuator so that it can be moved in and out 
of the sample chamber. In the Waiting for Sample state (5), the probe is fully 
inserted into the sample chamber so that the shaft of the probe is immersed in 

15 sample after the chamber is filled. The Measure Spectrum valve state (2) not only 
closes valves C and D, but also actuates both three-way valves controlling the 
linear actuator so that the probe is slowly extracted from the sample chamber 
while data is being collected. Upon completion of the Spectrum Collect 
operation, the sample in the sample chamber is ejected back into the sample 

20 transport line by opening valves C and E. 

Data Subsystem 

[0096] Referring again to Figure 2, the data subsystem includes an analyzer 302, 
which receives the response signal of the detector 106. The analyzer can be, for 

25 example, a computer capable of storing and processing the Raman data. Other 
functions of the analyzer can include, for example, developing the regression 
model and carrying out PCA/LWR analysis, as described below. In one 
embodiment described above, the data subsystem controls the motion of the 
sampling probe. In another embodiment described above, the data subsystem 

30 controls valves for filling and emptying the sample chamber. In another 
embodiment, the data subsystem compares the calculated value of one or more 
polymer properties to a target value, and adjusts one or more reactor parameters in 
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response to the deviation between calculated and target values. Reactor control is 
further described below. 

PCA/LWR Analysis 

5 [0097] The Raman spectrum includes information directly or indirectly related to 
various properties of the polyolefin sample. Conventionally, sample components 
are identified by the presence of unique spectral signatures, such as particular 
bands recognized as being due to particular vibrational modes of a molecule. 
Quantitative information such as concentration can then be obtained about a 
10 sample component by, for example, integrating the area under a particular peak 
and comparing the area to a calibration sample, by monitoring scattered intensity 
at a particular peak as a function of time, etc. In contrast to these conventional 
approaches, the present inventors have surprisingly found that polymer properties 
can be determined from Raman spectra without the need to identify or select 
15 particular spectral features, by using a multivariate model to correlate polymer 
properties with Raman scattering data. The model uses large, contiguous regions 
of the spectrum, rather than discrete spectral bands, thereby capturing large 
amounts of information density unavailable and unrecognized in conventional 
analysis. Further, the spectral data are correlated to polymer properties such as 
20 melt flow rates (defined below), densities, molecular weight distributions, etc., 
that are not readily apparent from optical spectra. 

[0098] In one embodiment, the data analysis described below is used to build and 
apply a predictive model for at least one property of the polyolefin particles 
selected from melt flow rate, density, molecular weight, molecular weight 

25 distribution, and functions thereof. 

[0099] As used herein, the term "melt flow rate" indicates any of the various 
quantities defined according to ASTM D-1238, including I 2 .i 6 , the melt flow rate 
of the polymer measured according to ASTM D-1238, condition E (2.16 kg load, 
190 °C), commonly termed the "melt index", and I 2 i. 6j the melt flow rate of the 

30 polymer measured according to ASTM D-1238, condition F (21.6 kg load, 
190 °C), commonly termed the "flow index." Other melt flow rates can be 
specified at different temperatures or different loads. The ratio of two melt flow 
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rates is the "Melt Flow Ratio" or MFR, and is most commonly the ratio of 
I21VI2.16. "MFR" can be used generally to indicate a ratio of melt flow rates 
measured at a higher load (numerator) to a lower load (denominator). 
[0100] As used herein, "molecular weight" indicates any of the moments of the 
5 molecular weight distribution, such as the number average, weight average, or Z- 
average molecular weights, and "molecular weight distribution" indicates the ratio 
of two such molecular weights. In general, molecular weights M can be computed 
from the expression: 

10 [0101] where Ni is the number of molecules having a molecular weight Mj. When 
n = 0, M is the number average molecular weight Mn. When n = 1, M is the 
weight average molecular weight Mw. When n = 2, M is the Z-average molecular 
weight Mz. These and higher moments are included in the term "molecular 
weight." The desired molecular weight distribution (MWD) function (such as, for 

15 example, Mw/Mn or Mz/Mw) is the ratio of the corresponding M values. 
Measurement of M and MWD by conventional methods such as gel permeation 
chromatography is well known in the art and is discussed in more detail in, for 
example, Slade, P. E. Ed., Polymer Molecular Weights Part II, Marcel Dekker, 
Inc., NY, (1975) 287-368; Rodriguez, F., Principles of Polymer Systems 3rd ed. y 

20 Hemisphere Pub. Corp., NY, (1989) 155-160; U.S. Patent No. 4,540,753; 
Verstrate et al. s Macromolecules, vol. 21, (1988) 3360; and references cited 
therein. 

[0102] Methods of the invention include obtaining a regression model for 
determining a polymer property, the regression model including principal 
25 component loadings and principal component scores; acquiring a Raman spectrum 
of a polyolefin sample; calculating a new principal component score from at least 
a portion of the Raman spectrum and the principal component loadings; and 
calculating the polymer property by applying the new principal component score 
to the regression model. 
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[0103] The regression model is preferable a locally weighted regression (LWR) 
model, using principal component analysis (PCA) eigenvectors. PCA is a well- 
known analytical method, and is described, for example, in Pirouette™ 
Multivariate Data Analysis for Windows software manual, Infometrix, Inc, 
5 Woodinville, WA (1985-2000), PLS_Toolbox™ software manual, Eigenvector 
Research, Inc., Manson, WA (1998), and references cited therein. LWR is 
described, for example, in Naes and Isaksson, Analytical Chemistry, 62, 664-673 
(1990), Sekulic et al., Analytical Chemistry, 65, 835A-845A (1993), and 
references cited therein. 
10 [0104] Principal Components Analysis is a mathematical method which forms 
linear combinations of raw variables to construct a set of mutually orthogonal 
eigenvectors (principal component loadings). Since the eigenvectors are mutually 
orthogonal, these new variables are uncorrelated. Further, PCA can calculate the 
eigenvectors in order of decreasing variance. Although the analysis computes a 
number of eigenvectors equal to the number of original variables, in practice, the 
first few eigenvectors capture a large amount of the sample variance. Thus, only a 
relatively small number of eigenvectors is needed to adequately capture the 
variance, and a large number of eigenvectors capturing minimal variance can be 
disregarded, if desired. 
20 [0105] The data are expressed in an m (row) by n (column) matrix X, with each 
sample being a row and each variable a column optionally mean centered, 
autoscaled, scaled by another function or not scaled. The covariance of the data 
matrix, cov(X), can be expressed as: 

cov(X) = X T X/(m-l) 

25 [0106] where the superscript T indicates the transpose matrix. The PCA analysis 
decomposes the data matrix as a linear combination of principal component scores 
vectors S,- and principal component loading vectors (eigenvectors) Lj, as follows: 

X = S,Li T + S 2 L2 T + S3L 2 T + ... 
[0107] The eigenvectors L; are eigenvectors of the covariance matrix, with the 

30 corresponding eigenvalues Xj indicating the relative amount of covariance 
captured by each eigenvector. Thus, the linear combination can be truncated after 
the sum of the remaining eigenvalues reaches an acceptably small value. 
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[0108] A model can be constructed correlating the Raman scattering intensity 
■with a polymer property in PCA space using various linear or nonlinear 
mathematical models, such as principal components regression (PCR), partial 
least squares (PLS), projection pursuit regression (PPR), alternating conditional 
5. expectations (ACE), multivariate adaptive regression splines (MARS), and neural 
networks (NN), to name a few. 

[0109] In a particular embodiment, the model is a locally weighted regression 
model. Locally Weighted Regression (LWR) assumes that a smooth non-linear 
function can be approximated by a linear or relatively simple non-linear (such as 
10 quadratic) function, with only the closest data points being used in the regression. 
The q closest points are used and are weighted by proximity, and the regression 
model is applied to the locally weighted values. 

[0110] In the calibration phase, Raman spectra are acquired, and the polymer 
properties of the sample are measured in the laboratory. The properties measured 

15 include those that the model will predict, such as density, melt flow rates, 
molecular weights, molecular weight distributions, and functions thereof. For a 
desired polymer property, the data set including the measured polymer properties 
the samples and the Raman spectral data for the samples is decomposed into PCA 
space to obtain a calibration data set. No particular number of calibration samples 

20 is required. One skilled in the art can determine the appropriate number of 
1 calibration samples based on the performance of the model and the incremental 
change in performance with additional calibration data. Similarly, there is no 
particular number of PCA eigenvectors required, and one skilled in the art can 
choose an appropriate number based on the amount of variance captured a 

25 selected number of eigenvectors and the incremental effect of additional 
eigenvectors. 

[0111] The LWR model can be validated using methods known in the art. It is 
convenient to divide the calibration samples into two sets: a calibration data set, 
and a validation data set. The calibration data set is used to develop the model, 
30 and to predict the appropriate polymer property for the samples in the validation 
data set, using the validation data set Raman spectra. Since the chosen polymer 
property for the validation data set samples is both calculated and measured, the 
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effectiveness of the model can be evaluated by comparing the calculated and 
measured values. 

[0112] The validated model can then be applied to sample spectra to predict the 
desired polymer property or properties. 
5 [01 13 J If desired, a single model can be used to predict two or more polymer 
properties. Preferably, separate models are developed for e&ch polymer property. 
Thus, in one embodiment, the present invention includes: obtaining a first 
regression model for determining a first polymer property, the first regression 
model including first principal component loadings and first principal component 

10 scores; obtaining a second regression model for determining a second polymer 
property, the second regression model including second principal component 
loadings and second principal component scores; acquiring a Raman spectrum of 
a sample comprising polyolefin; calculating a new first principal component score 
from at least a portion of the Raman spectrum and the first principal component 

15 loadings; calculating a new second principal component score from at least a 
portion of the Raman spectrum and the second principal component loadings; 
calculating the first polymer property by applying the new first principal 
component score to the first regression model; and calculating the second polymer 
property by applying the new second principal component score to the second 

20 regression model. 

[0114] Of course, more than two polymer properties can be determined by 
including third or more regression models. Advantageously, multiple polymer 
properties can be determined essentially simultaneously by using the same Raman 
spectrum and applying several regression models to the spectral data. 

25 [0115] In a particular embodiment, two regression models are used, and both a 
melt flow rate (such as melt index h.ie or flow index I21.6) and density are 
determined. 

Reaction Control 

30 [0116] In one embodiment, the calculated polymer property is compared to a 
target polymer property, and at least one reactor parameter is adjusted based on 
the deviation between the calculated and target polymer property. The at least one 
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reactor parameter can include the amounts of monomer, comonomer, catalyst and 
cocatalyst, the operating temperature of the reactor, the ratio of comonomer(s) to 
monomer, the ratio of hydrogen to monomer or comonomer, and other parameters 
that affect the chosen polymer property. For example, if the chosen polymer 
5 property is density and the density calculated from the PCA/LWR model is lower 
than a target density, a reactor parameter can be adjusted to increase density, such 
as, for example, reducing the comonomer feed rate and/or increasing the monomer 
feed rate. 

[0117] For example, in the case of the fluidized bed polymerization of olefins, 
10 hydrogen can serve as a chain transfer agent. In this way, the molecular weight of 
the polymer product can be controlled. Additionally, varying the hydrogen 
concentration in olefin polymerization reactors can also vary the polymer melt 
flow rate, such as the melt index I 2 .i 6 (MI). The present invention allows control of 
the reactor to produce polymer having a selected MI range. This is accomplished 
15 by knowing the relationship between hydrogen concentration and the MI of 
polymers produced by a specific reactor, and programming the target MI or MI 
range into a reactor control system processor. By monitoring the polymer MI data 
generated by the Raman analyzer and comparing this data to the target MI range, 
the flow of hydrogen into the reactor vessel may be adjusted so that the MI range 
20 of the polymer product remains within the target MI range. 

[0118] It will be understood by those skilled in the art that other reactor 
constituent properties and other reactor parameters can be used. In a similar way 
as described above, the final polymer properties may be achieved by controlled 
metering reactor parameters in response to data generated by the Raman analyzer. 

25 

EXAMPLES 

[0119] Laboratory determinations of density (g/cm 3 ) used a compression molded 
sample, cooled at 15 °C per hour and conditioned for 40 hours at room 
temperature according to ASTM D1505 and ASTM D1928, procedure C. 
30 [0120] Laboratory determinations of melt flow rates were carried out at 190 °C 
according to ASTM D-1238. I 2 i.6 is the "flow index" or melt flow rate of the 
polymer measured according to ASTM D-1238, condition F, and I 2 .i6 is the "melt 
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index" or melt flow rate of the polymer measured according to ASTM D-1238, 
condition E. The ratio of I 2 i. 6 to I Zi6 is the "melt flow ratio" or "MFR". 
[0121] EXCEED™ 350 is a gas-phase metallocene produced LLDPE 
ethylene/hexene copolymer with a Melt Index (Izifi) of 1.0 g/10 min, and a density 
5 of 0.918 g/cm 3 , available from ExxonMobil Chemical Co., Houston, TX. The 
EXCEED™ 350 resin is now marketed as EXCEED™ 3518. 
[0122] EXCEED™ 357 is a gas-phase metallocene produced LLDPE 
ethylene/hexene copolymer with a Melt Index (I 2 . I6 ) of 3.4 g/10 min, and a density 
of 0.917 g/cm 3 , available from ExxonMobil Chemical Co., Houston, TX. The 
10 EXCEED™ 357 resin is now marketed as EXCEED™ 35 1 8. 

[0123] ExxonMobil LL-1002 is a gas-phase Ziegler-Natta produced LLDPE 
ethylene/butene copolymer resin having a Melt Index (I 2 .i 6 ) of 2.0 g/10 min, and a 
density of 0.918 g/cm 3 , available from ExxonMobil Chemical Co., Houston, TX. 
[0124] ExxonMobil LL-1 107 is a gas-phase Ziegler-Natta produced LLDPE 
1 5 ethylene/butene copolymer resin having a Melt Index (I 2 .i 6 ) of 0.8 g/1 0 min, and a 
density of 0.922 g/cm 3 , available from ExxonMobil Chemical Co., Houston, TX. 
[0125] ExxonMobil LL-6100 is a gas-phase Ziegler-Natta produced LLDPE 
ethylene/butene copolymer resin having a Melt Index (I 2 ., 6 ) of 20 g/10 min, and a 
density of 0.925 g/cm 3 , available from ExxonMobil Chemical Co., Houston, TX. 
20 [0126] ExxonMobil LL-6101 is a gas-phase Ziegler-Natta produced LLDPE 
ethylene/butene copolymer resin having a Melt Index (I 2 ., 6 ) of 20 g/10 min, and a 
density of 0.925 g/cm 3 , available from ExxonMobil Chemical Co., Houston, TX. 
[0127] ExxonMobil LL-6201 is a gas-phase Ziegler-Natta produced LLDPE 
ethylene/butene copolymer resin having a Melt Index (I 2 .i 6 ) of 50 g/10 min, and a 
25 density of 0.926 g/cm 3 , available from ExxonMobil Chemical Co., Houston, TX. 

Examples 1-3 

[0128] Examples 1-3 were used to show the feasibility of embodiments of the 
invention. In Examples 1-3, measurements were made in the laboratory, 
30 simulating the measurements that would be made on-line in a polymerization 
reactor. 
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[0129] The Raman system used for Examples 1-3 was a Kaiser Optical Holoprobe 
Process Raman Analyzer, available from Kaiser Optical Systems, Inc., Ann 
Arbor, Michigan. The Raman system used a 125 mW diode laser operating at 785 
nm, and was equipped with a probe with 2.5 (6.3 cm) inch imaging optics fiber- 
5 optically coupled to the instrument, a holographic notch filter, holographic 
dispersion grating, cooled CCD detector (-40 °C), and computer for analyzer 
control and data analysis. A more complete description of this commercial 
instrument can be found in "Electro-Optic, Integrated Optic, and Electronic 
Technologies for Online Chemical Process Monitoring," Proceedings SPIE, vol. 
10 3537, pp. 200-212 (1998), the disclosure of which is incorporated herein by 
reference for purposes of U.S. patent practice. 

[0130] Data collection was accomplished by positioning the Raman probe above 
the surface of a polymer granule sample at a distance of about 2.5 inches (6.3 cm). 
The probe was fiber optically coupled to the Raman analyzer for both excitation 

15 and scattering signals. Data were collected from each sample for three minutes 
(i.e., signal averaged for 3 minutes). The CCD detector is sensitive to cosmic 
rays, which can cause spurious signals in array elements. "Cosmic ray checking" 
is a detector function that checks for these artifacts and discards them. In the 
following examples, the cosmic ray checking function was used. 

20 [0131] Raman spectra were collected over the region of 100 to 3500 cm" 1 . Three 
consecutive spectra were collected for each sample used. The samples were 
obtained from either of two gas-phase fluidized bed reactors producing 
copolymers of ethylene and butene or hexene, using metallocene catalysts. 
Laboratory measurements of melt index and/or density were also made for each 

25 sample. 

[0132] The data were divided into calibration sets, used to develop the PCA/LWR 
models, and validation sets, used to evaluate the accuracy of the model. Separate 
models were developed for a relatively low melt index range, a relatively high 
melt index range, and density. 
30 Example 1: Low Melt Index Model 

[0133] Seventy-three polymer samples were evaluated. The samples were divided 
into a group of 50 used for calibration (model development) and a group of 23 
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used for model validation. Each sample was a metallocene-catalyzed LLDPE 
resin, with hexene comonomer, in a melt index range of from about 0.6 to about 
1.2 g/10 min. Raman spectra and laboratory melt index measurements were 
collected as described above. 
5 [0134] The lab values of melt index and the Raman spectra of the calibration data 
set were used to create a locally-weighted regression model for low range melt 
index, using principal component loadings and principal component scores. The 
measured melt indexes, predicted melt indexes, and deviations (i.e., deviation of 
the actual melt index from the prediction of the LWR model) are shown in Table 
10 1. 
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Table 1 : Low MI Calibration 



MI (Lab) 
(dg/min) 


MI (Model) 
(dg/min) 


AMI (a) 
(dg/min) 


MI (Lab) 
(dg/min) 


MI (Model) 
(dg/min) 


AMI< B) 
(dg/min) 


0.678 


0.663401 


-0.0146 


0.82 


0.823629 


0.003629 


0.678 


0.675728 


-0.00227 


0.831 


0.841478 


0.010478 


0.679 


0.685591 


0.006591 


0.831 


0.8089 


-0.0221 


0.679 


0.653462 


-0.02554 


0.84 


0.819804 


-0.0202 


0.687 


0.699942 


0.012942 


0.84 


0.838078 


-0.00192 


0.687 


0.709433 


0.022433 


0.92 


0.934314 


0.014314 


0.696 


0.700481 


0.004481 


0.92 


0.93859 


0.01859 


0.696 


0.696309 


0.000309 


1.06 


1.049136 


-0.01086 


0.7 


0.689811 


-0.01019 


1.06 


1.07161 


0.01161 


0.7 


0.694658 


-0.00534 


1.07 


1.080271 


0.010271 


0.705 


0.690562 


-0.01444 


1.07 


1.079701 


0.009701 


0.705 


0.706591 


0.001591 


1.08 


1.090437 


0.010437 


0.706 


0.69476 


-0.01124 


1.08 


1.055101 


-0.0249 


0.706 f 


0.718346 


0.012346 


1.098 


1.117367 


0.019367 


0.714 


0.706535 


-0.00746 


1.098 


1.092972 


-0.00503 


0.714 


0.703178 


-0.01082 


1.1 


1.083835 


-0.01617 


0.7546 


0.786602 


0.032002 


1.1 


1.071211 


-0.02879 


U./DHO 


U. / /40I0 


0.020016 


1.11 


1.115756 


0.005756 


0.772 


0.781622 


0.009622 


1.11 


1.106827 


-0.00317 


0.772 


0.779611 


0.007611 


1.11 


1.085486 


-0.02451 


0.773 


0.775132 


0.002132 


1.11 


1.096664 


-0.01334 


0.773 


0.777378 


0.004378 


1.15 


1.142874 


-0.00713 


0.808 


0.800435 


-0.00757 


1.15 


1.1283 


-0.0217 


0.808 


0.824823 


0.016823 


1.1811 


1.200165 


0.019065 


0.82 


0.825021 


0.005021 


1.1811 


1.198869 


0.017769 



Model (predicted) MI minus Lab (measured) MI 



[0135] The Raman spectra of the validation data set were collected, and new 
5 principal component scores were calculated from the validation spectra. Using the 
locally-weighted regression model, the melt index of each validation sample was 
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then calculated. The measured melt indexes, predicted melt indexes, and 
deviations (i.e., deviation of the actual melt index from the prediction of the LWR 
model) are shown in Table 2. 



Table 2: Low MI Validation 



MI (Lab) 
(dg/min) 


MI (Model) 
(dg/min) 


AMI (a) 
(dg/min) 


MI (Lab) 
(dg/min) 


(dg/min) 


(dg/min) 


0.55 


0.561835 


0.011835 


0.883 


0.802692 


-0.08031 


0.55 


0.579349 


0.029349 


0.883 


0.856272 


-0.02673 


0.55 


0.57315 


0.02315 


0.883 


0.849839 


-0.03316 


0.616 


,0.654254 


0.038254 


0.95 


1.083123 


0.133123 


0.616 


0.637083 


0.021083 


0.95 


1.022883 


0.072883 


0.616 


0.667328 


0.051328 


0.95 


1.021329 


0.071329 


0.622 


0.6504 


0.0284 


1.065 


1.006358 


-0.05864 


0.622 


0.635863 


0.013863 


1.065 


0.950208 


-0.11479 


0.622 


0.669156 


0.047156 


1.065 


0.978949 


-0.08605 


0.679 


0.644011 


-0.03499 


1.142 


1.14752 


0.00552 


0.679 


0.632626 


-0.04637 


1.142 


1.12363 


-0.01837 


0.679 


0.634522 


-0.04448 









5 Model (predicted) MI minus Lab (measured) MI 



[0136] Figure 6A depicts the data from Tables 1 and 2 graphically. The line in 
the Figure is the model prediction. The calculated R 2 value was 0.99 for the 
calibration set, with a standard error of 0.0155, and 0.92 for the validation set, 
10 with a standard error of 0.059. 

Example 2: High Melt Index Model 

[0137] An analysis was carried out as in Example 1, using higher melt index 
samples. Thirty-four polymer samples were evaluated. These samples were used 
15 as calibration samples for model development, but a validation subset was not 
used. Each sample was a metallocene-catalyzed LLDPE resin, with butene 
comonomer, in a melt index range of from about 4 to about 60 g/10 min. Raman 
spectra and laboratory melt index measurements were collected as described 
above. 
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[0138] The lab values of melt index and the Raman spectra of the calibration data 
set were used to create a locaUy-weighted regression model for high range melt 
index, using principal component loadings and principal component scores. The 
measured melt indexes, predicted melt indexes, and deviations (i.e., deviation of 
5 the actual melt index from the prediction of the LWR model) are shown in Table 
3. 



Table 3: High MI Calibration 



MI (Lab) 
(dg/min) 


MI (Model) 
(dg/min) 


AMI (a) 
(dg/min) 


MI (Lab) 
(dg/min) 


MI flVlodeH 
(dg/min) 


AMT (o) 

LALVll 

(dg/min) 


4.341 


4.513 


0.172 


30.68 


29.118 


-1.562 


4.341 


4.467 


0.126 


32.93 


32.112 


-0.818 


8.613 


8.433 


-0.18 


32.93 


32.459 


-0.471 


8.613 


8.314 


-0.299 


33.68 


34.658 


0.978 


10.499 


9.978 


-0.521 


33.68 


34.233 


0.553 


10.499 


10.768 


0.269 


36.6 


37.216 


0.616 


12.547 


13.013 


0.466 


36.6 


36.989 


0.389 


12.547 


12.971 


0.424 


45.15 


44.433 


-0.717 


18.61 


17.955 


-0.655 


45.15 


45.001 


-0.149 


18.61 


17.885 


-0.725 


48.07 


48.966 


0.896 


19.81 


21.009 


1.199 


48.07 


49.207 


1.137 


19.81 


20.893 


1.083 


51.41 


49.879 


-1.531 


21.59 


22.011 


0.421 


51.41 


50.554 


-0.856 


21.59 


22.314 


0.724 


55.56 


57.213 


1.653 


22.79 


22.291 


-0.499 


55.56 


56.667 


1.107 


22.79 


23.109 


0.319 


57.41 


57.942 


0.532 


30.68 


28.212 


-2.468 


57.41 


58.217 


0.807 



Model (predicted) MI minus Lab (measured) MI 



10 [0139] Figure 6B depicts the data from Table 3 graphically. The line in the 
Figure is the model prediction. The calculated R 2 value was 0.99, with a standard 
error of 0.91. 
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Example 3: Density Model 

[0140] An analysis was carried out as in Example 1, using density rather than 
melt index as the predicted property. A subset of 22 of the polymer samples used 
in Example 1 were evaluated. These samples were used as calibration samples for 
5 model development, but a validation subset was not used. Each sample was a 
metallocene-catalyzed LLDPE resin, with hexene comonomer. Raman spectra 
and laboratory density measurements were collected as described above. 
[0141] The lab values of density and the Raman spectra of the calibration data set 
were used to create a locally-weighted regression model for density, using 
10 principal component loadings and principal component scores. The measured 
densities, predicted densities, and deviations (i.e., deviation of the actual density 
from the prediction of the LWR model) are shown in Table 4. 



Table 4: Density Calibration 



p(Lab) 
(g/cm 3 ) 


p (Model) 
(g/cm 3 ) 


(g/cm 3 ) 


p(Lab) 
(g/cm 3 ) 


p (Model) 
(g/cm 3 ) 


P(8> 3 

(g/cm 3 ) 


0.9183 


0.919018 


0.000718 


0.9218 


0.922121 


0.000321 


. 0.9183 


0.919053 


0.000753 


0.922 


0.921901 


-0.000099 


0.9185 


0.917859 


-0.00064 


0.922 


0.922797 


0.000797 


0.9185 


0.917786 


-0.00071 


0.9226 


0.921872 


-0.00073 


0.9195 


0.918575 


-0.00092 


0.9226 


0.922369 


-0.00023 


0.9195 


0.918499 


-0.001 


0.9244 


0.924316 


-0.000084 


0.9196 


0.919342 


-0.00026 


0.9244 


0.924075 


-0.00033 


0.9196 


0.919943 


0.000343 


0.9249 


0.924893 


-0.000007 


0.9212 


0.921674 


0.000474 


0.9249 


0.924031 


-0.00087 


0.9212 


0.921701 


0.000501 


0.9262 


0.926252 


0.000052 


0.9218 


0.92193 


0.00013 


0.9262 


0.925936 


-0.00026 



Model (predicted) density minus Lab (measured) density 



15 

[0142] Figure 7 depicts the data graphically. The line in the Figure is the model 
prediction. The calculated R 2 value was 0.95, with a standard error of 0.00057. 
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Examples 4-5 

[0143] Examples 4-5 demonstrate the effectiveness of the inventive methods on- 
line in a polymerization reaction system, for melt index determination. 
5 [0144] The Raman system used for Examples 4-5 was as described for Examples 
1-3, except that the laser was a 200 mW mode-stabilized diode laser operating at 
785 nm. Polymer samples from either of two gas-phase fluidized-bed reactors 
were taken using the sampling system described above. 

[0145] The data were divided into calibration sets, used to develop the PCA/LWR 
10 models, and validation sets, used to evaluate the accuracy of the model. Separate 
models were developed for a melt index (Examples 4-5) and density (Examples 6- 
7). In addition, separate models were developed for each of the two gas-phase 
reactors. The two reactors are denoted "Reactor 1" and "Reactor 2" below. 

15 Example 4: Melt Index Model Reactor 1 

[0146] Two hundred eighty-five polymer samples were evaluated. The samples 
were divided into a group of 216 used for calibration (model development) and a 
group of 69 used for model validation. Each sample was a metallocene-catalyzed 
LLDPE resin, in a melt index range of from less than 1 to about 15 g/10 min. 

20 Raman spectra and laboratory melt index measurements were collected as 
described above. 

[0147] The lab values of melt index and the Raman spectra of the calibration data 
set were used to create a locally-weighted regression model for melt index, using 
principal component loadings and principal component scores. The measured 

25 melt indexes and predicted melt indexes are shown in Tables 5A-5B. The 
deviations are not shown in the table, but are readily calculated from the tabulated 
data. The data are shown in the order taken (by column, within each table), to 
illustrate the effectiveness of the model under changing polymer conditions. A 
symbol "Vn" before an entry indicates that the nth set of validation spectra were 

30 taken before the marked entry, as shown by the corresponding notation in Table 6. 
Table 5B is a continuation of Table 5 A. 
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Table 5A: MI Calibration, Reactor 1 



M7 CI rM 
(dg/min) 


mj (lYLoueij 
(dg/min) 


Ml (Lab; 
(dg/min) 


Ml (Model) 
(dg/min) 


MI (Lab) 

( d&Anin^ 
mill j 


MI (Model) 

Inn /tyi mi 


4.997 


5.013 


3.506 


3.440 




0.973 


4.413 


4.390 


3.554 


3.401 




0.961 


4.559 


4.410 


3.554 


3.474 


W.7JO 


U.y I / 


3.511 


3.633 


3.541 


3.540 




0.940 


3.481 


3.521 


3.576 


3.713 


ft 073 




3.315 


3.391 


3.576 


3.679 


ft QA& 




3.301 


3.286 


3.630 


3.664 


ft Q70 


u.you 


3.369 


3.211 


3.630 


3.664 


1 1 ^ 
1 . 1 JJ 


l.UoJ 


3.460 


3.607 


3.626 


3.563 


1 .loo 




3.391 


3.481 


3.618 


3.652 


1 1 20 
LIB/ 


1.231 


3.380 


3.301 


3.346 


3.257 




1. 1U4 


3.523 


3.629 


3.409 


3 399 


1 1 38 

l.lJO 


1.193 


3.370 


3.294 


3.409 


3.426 




u.yyo 


3.537 


3.522 


(VI> 3.411 


3.342 


ft Q77 


u.yo/ 


3.534 


3.559 


3.572 


3.743 




u.y /J 


3.432 


3.407 


2.351 


2.402 


0 Q7ft 




3.518 


3.671 


1.544 


1.574 


ft QRS 




3.555 


3.562 


1.348 


1.364 


ft QS? 


V.yOZ, 


3.380 


3.299 


1.163 


1.140 




n 01 c 
\J.y l o 


3.320 


3.308 


1.106 


1.095 


0 071 


u.yuu 


3.470 


3.523 


1.072 


1.100 


1.017 


v.yo i 


3.380 


3.405 


1.098 


1.103 


1 t\S\J*/ 


1 ft£1 


3.380 


3.277 


1.071 


1.110 


1 010 


i ni7 

1 .U 1Z 


3.370 


3.328 


0.987 


0.971 


1.030 


1 07S 


3.370 


3.400 


1.009 


0.994 


0.986 


n 070 

\j.y /y 


3.354 


3.290 


1.005 


0.998 


0.944 


ft 017 


3.354 


3.540 


0.978 


0.980 


0.947 


0.953 


3.523 


3.327 


1.009 


1.010 


0.955 


0.972 


3.523 


3.473 


0.991 


1.039 


0.932 


0.928 


3.491 


3.551 


1.002 


0.991 


0.947 


0.944 


3.582 


3.613 


1.038 


1.094 


0.984 


0.990 


3.582 


3.612 


1.035 


1.000 


0.972 


0.994 


3.493 


3.464 


1.016 


1.023 


0.998 


0.990 


3.493 


3.375 


0.940 


0.932 


0.991 


1.039 


3.523 


3.596 


0.970 


0.980 


1.060 


1.004 


3.506 


3.483 


0.980 


0.979 


1.041 


1.013 
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Table 5B: MI Calibration, Reactor 1 , continued 



Ml (Lab) 
(dg/min) 


MI (Model! 
(dg/min) 


MT CI ah\ 
(dg/min) 


Ml (Model » 
(dg/min) 


TV/fT fx n u,\ 

(dg/min) 


Mi (JVioaeij 
(dg/min) 


0.989 


1.000 


0.938 


0.942 


3.490 


3.311 


0.921 


0.910 


0.988 


0.960 


3.528 


3.466 


0.908 


0.880 


1.006 


1.039 


3.538 


3.555 


0.951 


0.962 


0.982 


1.002 


3.592 


3.403 


0.976 


0.990 


0.946 


0.978 


3.372 


3.441 


0.965 


0.941 


0.964 


0.930 


w 3.580 


3.726 


0.970 


0.992 


1.010 


0.962 


3.259 


3.109 


0.966 


1.002 


1.030 


1.082 


3.302 


3.319 


0.998 


1.092 


1.040 


1.019 


3.437 


3.572 


0.983 


0.963 


1.080 


1.103 


3.397 


3.429 


0.985 


0.973 


1.020 


1.031 


3.449 


3.401 


0.990 


0.999 


1.040 


1.039 


3.513 


3.329 


0.990 


1.024 


™ 1.061 


1.092 


3.771 


3.702 


0.993 


1.003 


1.546 


1.552 


3.986 


3.918 


0.968 


0.982 


2.043 


1.993 


4.575 


4.428 


0.997 


0.971 _ 


2.381 


2.402 


5.000 


4.892 


1.006 


0.982 


2.751 


2.772 


6.054 


6.203 


^ 0.939 


0.926 


3.054 


2.994 


7.452 


7.624 


0.966 


0.953 


3.414 


3.540 


8.392 


8.012 


0.992 


1.017 


3.342 


3.254 


10.630 


10.171 


0.989 


0.993 


3.550 


3.453 


10.630 


10.171 


0.951 


0.937 


3.580 


3.429 


12.530 


12.779 


i.030 


1.018 


3.550 


3.445 


13.110 


13.671 


1.000 . 


1.005 


3.610 


3.454 


13.879 


13.698 


0.959 


0.939 


3.528 


3.310 


13.952 


13.498 


0.954 


0.957 


3.24'6 


3.152 


13.627 


13.593 


0.940 


0.910 


3.523 


3.391 


13.295 


12.998 


0.985 


0.991 


3.620 


3.662 


13.393 


. 13.876 


0.980 


0.991 


3.691 


3.604 


13.146 


13.029 


0.955 


0.931 


3.713 


3.700 


12.810 


13.014 


0.930 


0.909 


3.451 


3.619 


11.989 


11.903 


0.910 


0.891 


3.439 


3.293 


10.670 


11.003 


0.940 


0.963 


3.501 


3.701 


12.181 


12.292 


0.980 


1.003 


3.263 


3.331 


12.711 


12.625 


0.980 


0.993 


3.433 


3.383 


1.120 


1.204 


0.960 


0.968 


3.477 


3.579 


1.002 


1.002 
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[0148] The Raman spectra of the validation data set were also collected, and new 
principal component scores were calculated from the validation spectra. Using the 
locally-weighted regression model, the melt index of each validation sample was 
then calculated. The measured and predicted melt indexes are shown in Table 6. 
5 Acquisition of the validation spectra was interspersed with acquisition of the 
calibration spectra, at the corresponding "Vn" positions. 



Table 6: MI Validation, Reactor 1 



MI (Lab) 
(dg/min) 


MI (Model) 
(dg/min) 


Ml (Lab) 
(dg/min) 


MI (Model) 
(dg/min) 


MI (Lab) 
(dg/min) 


MI (Model) 
(dg/min) 


VI: 3.471 


. 3.514 


0.965 


0.991 


0.933 


0.960 


3.443 


3.503 


1.000 


0.973 


0.859 


0.811 


3.438 


3.371 


0.995 


0.999 


0.934 


0.903 


3.493 


3.421 . 


0.964 


0.952 


0.980 


1.011 


3.417 


3.561 


0.934 


0.943 


0.910 


0.880 


3.354 


3.365 


0.946 


0.967 


0.890 


0.920 


3.454 


3.604 


0.943 


0.92O 


0.900 


0.899 


3.531 


3.594 


0.928 


0.892 


0.970 


0.992 


3.557 


3.500 


0.931 


0.950 


0.980 


0.962 


3.521 


3.498 


0.967 


0.972 


0.990 


1.053 


3.440 


3.352 


0.949 


0.957 


V4: 3.470 


3.625 


3.507 


3.521 


1.025 


0.980 


3.620 


3.772 


3.596 


3.569 


1.029 


1.089 


3.420 


3.400 


3.659 


3.623 


1.032 


1.012 


3.504 


3.387 


3.554 


3.648 


1.025 


1.034 


3.682 


3.598 


3.565 


3.662 


0.999 


1.004 


3.597 


3.784 


3.605 


3.807 


0.995 


0.970 


3.531 


3.724 


3.573 


3.531 


1.009 


0.998 


3.399 


3.412 


3.456 


3.604 


1.035 


1.029 


3.590 


3.498 


3.501 


3.586 


1.048 


1.011 


3.520 


3.500 


3.500 


3.398 


1.012 


1.029 


3.431 


3.548 


V2: 0.980 


0.998 


V3: 1.095 


1.118 


3.391 


3.293 


0.966 


0.982 


1.114 


1.092 


3.288 


3.412 
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[0149] Figure 8A depicts the data from Tables 5A, 5B and 6 graphically. The line 
in the Figure is the model prediction. The calculated R 2 value was 0.999, with a 
standard error of 2.78%. 

5 Example 5: Melt Index Model. Reactor 2 

[0150] The procedure described in Example 4 was followed, except as noted, 
sampling this time from the Reactor 2 polymer. Two hundred ninety-one polymer 
samples were evaluated. The samples were divided into a group of 266 used for 
calibration (model development) and a group of 25 used for model validation. 
10 Each sample was a Ziegler-Natta-catalyzed LLDPE resin, in a melt index range of 
from less than 1 to about 60 g/10 min. Raman spectra and laboratory melt index 
measurements were collected as described above. 

[0151] The lab values of melt index and the Raman spectra of the calibration data 
set were used to create a locally-weighted regression model for melt index, using 

15 principal component loadings and principal component scores. The measured 
melt indexes and predicted melt indexes are shown in Tables 7A-7B. The 
deviations are not shown in the table, but are readily calculated from the tabulated 
data. The data are shown in the order taken (by column, within each table), to 
illustrate the effectiveness of the model under changing polymer conditions. A 

20 symbol "Vn" before an entry indicates that the nth set of validation spectra were 
taken before the marked entry, as shown by the corresponding notation in Table 8. 
Table 7B is a continuation of Table 7A. In Tables 7A and 7B, the units of melt 
index (MI) are dg/min. 
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Table 7A: MI Calibration, Reactor 2 



MI 

(Lab) 


tun 
(Model) 


MT 

MA 

(Lab) 


Ml 

(Model) 


IV if! 

MI 
(Lab) 


MI 
(Model) 


MI 


MI 


| 0,678 


0.669 


0.890 


0.862 


53.370 




91 T CA 


21.888 


2.008 


2.105 


0.837 


0.843 


52.750 


52.559 


91 9** A 


on ii<:<r 
20.405 


1.410 


1.376 


V1 17.580 


18.001 


50.660 


51.032 






0.992 


0.988 


19.194 


18.899 


51.720 


50.759 


9ft 1*7 ft 


On vlO**7 

ZU.4/7 


0.832 


0.859 


20.280 


21.021 


48.530 


47.667 


iy,3r ft 


zu.u / / 


0.758 


0.780 


19.588 


20.091 


44.160 


45.221 


22 0*4 n 


99 V7A 


0.712 


0.688 


19.732 


18.979 


47.000 


46.821 


1 .ZJU 




0.673 


0.690 


20.910 


22.098 


53.370 


54.202 


*#- f 7u 


*3 flfM 
J.OUl 


0.670 


0.710 


20.070 


20.473 


41.750 


40.512 






0.721 


0.690 


19.800 


19.070 


48.360 


49.848 




1 ooc 
1.770 


0.753 


0.774 


20.900 


20.078 


50.890 


49. 1 11 


9 2QB 




0.751 


0.779 


22.080 


21.874 


43.810 


43.084 


1 047 


1 J25J9 
A. OO / 


0.780 


0.810 


20.080 


19.659 


43.850 


44.106 




1 fi£7 


0.811 


0.779 


19.616 


19,223 


46.200 


47 485 


9 fKQ 




0.792 


0.81O 


19.829 


19.629 


48.220 


48 944 


9 1^1 




0.782 


0.750 


37.090 


17.651 


49.950 


49 004 




4.H7 


0.753 


0.775 


18.086 


17.888 


49.590 


50.672 


9 17H 


1 AG^ 


0.767 


0.798 


17.638 


16.844 


^41.540 


41.094 


1 A on 


9 1 77 


0.706 


0.700 


18.637 


17.974 


21.320 


22.445 




9 914^ 


0.878 


0.892 


20.010 


20.119 


17.983 


18.241 


2 1 10 


i 


0,858 


0.823 


19,568 


19.103 


17.233 


16.869 


2 A*ifi. 


1,707 


0.817 


0.829 


V2 27.270 


27.906 


19.677 


19-311 


2,000 


1 GOO 


0.857 


0.802 


42.780 


43.099 


19.063 


19.921 


1 Q17 

1.71 / 


i ov4 ! 
1 i>V*» 


0.849 


0.836 


48.560 


47.956 


19,919 


19.107 


1 07A 


9 lAl 


0.779 


0.750 


51.950 


52.302 


21.510 


21.444 




9 AA1 


0.765 


0.770 


51.270 


51.037 


20.840 


20.77) 


2.077 


1.985 


0.742 


0.727 


49.950 


50.119 


20.500 


19.295 


2.035 


2.103 


0.806 


0.800 


45.610 


46.117 


20.230 


21.011 


2.007 


2.U0 


0.810 


0.801 


47.440 


46.938 


21.230 


20.659 


1.980 


2.004 


0.827 


0.839 


53.620 


52.476 


21.670 


20.997 


1.950 


1.891 


0.778 


0.789 


55.210 


54.998 


20.590 


21.264 


1.880 


1.871 


0.796 


0.763 


50.010 


49.483 


23.140 


22.7B4 


1,990 


2.109 


0.768 


0.712 


44.040 


44.884 


22.460 


21.997 


2.230 


2.190 


0.899 


0.912 


42.780 


41.009 


20.640 


20.883 


2.100 


1.962 


0.946 


0.963 


47.940 


47.444 


21.260 


20.799 


1.998 


2.119 


0.965 


1.002 


53.720 


52.798 


19.856 


20.231 


1.910 


1.967 
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Tabic 7B: MI Calibration, Reactor 2» continued 



Ml 
(Lab) 


MI 

(Model) 


MI 
(Lab) 


MI 
(Model) 


MI 

(Lab) 


MI 

(Model) 


u 

JVLl 

(Lab) 


Ml 

(Model) 


2.204 


2.187 


0.974 


0.980 


2.010 


1.992 


1.242 


1.219 


Z350 


2.410 


0.992 


1.006 


1.346 


1.375 


L320 


1.331 


Z201 


2.177 


0.924 


0.898 


0.945 


0.972 


1.396 


1387 


2.050 


1.939 


0.983 


0.978 


0.700 


0.700 


1.480 


1.520 


2.120 


2.098 


0-970 


0.952 


0.825 


0.830 


1.594 


1.554 


2.000 


2.079 


1.079 


1.108 


0.843 


0,852 


1.525 


1.501 


2.040 


2.101 


1.077 


0.997 


0.792 


0.804 


1.576 


1.629 


2.100 


2.008 


1.093 


1.121 


0.791 


0.801 


1.664 


1.711 


2.040 


2.U3 


1.108 


1.110 


0.796 


0.799 


1.544 


1.557 


1.950 


1.889 


1.071 


0.997 


0.745 


0.756 


1.962 


1.891 


V5 0.945 


0.902 


1.013 


0.999 


0.777 


0.791 


5.093 


4.894 


0970 


0.978 


0.980 


0.980 


0.734 


0.720 


9.130 


9-297 


0.965 


0.954 


1.061 


1.046 


0.711 


0.720 


12.063 


U.999 


1.281 


1.299 


1.005 


1.018 


0.763 


0.777 


13.815 


13.684 


1.445 


1.455 


0.961 


0.967 


0.778 


0,801 


13.262 


13.555 


1,502 


1.552 


1.005 


0.997 


0685 


0.673 


16.134 


16.643 


1.373 


1.299 


0.980 


0.977 


0.769 


0.776 


15.180 


15322 


1.365 


1399 ! 


0.864 


0.877 


0-760 


0.743 


15.845 


16,015 


1.420 


1.390 


0.891 


0-903 


0.738 


0,750 


15.730 


15.926 


1.462 


1.442 


0.996 


1.043 


0.726 


0.701 


12.221 


12.442 


1.674 


1,739 


1.054 


1.044 


0.719 


0.742 


11.531 


11.735 


1.863 


1.920 


1.017 


1.009 


0.706 


0.688 


12.532 


12.221 


2.168 


2.122 


1.023 


0.995 


0.781 


0.743 


12.358 


12.471 


1.979 


1.948 


2.079 


1.997 


0.797 


0.822 


12.538 


12.882 


3.279 


3.309 


1.963 


2,047 


0.750 


0.770 


12,549 


12.555 


0.969 


1.002 


1.963 


2.065 


0.779 


0.788 


12.948 


12.507 


1.018 


1.040 


1.880 


1.841 


0.806 


0.810 


13.413 


13.119 


1.078 


1.039 


2.070 


2.109 


0.768 


0.773 


12.543 


12.629 


1.009 


0.988 


2.180 


2.116 


0.896 


0.910 


12.500 


12.409 


1.034 


0.992 


2.290 


2.341 


1.142 


1.172 


12.111 


12.427 


1.005 


1.056 


2.150 


2.098 


1.176 


1.149 


11.957 


11.883 1 



[0152] The Raman spectra of the validation data set were also collected, and new 
principal component scores were calculated from the validation spectra. Using the 
5 locally-weighted regression model, the melt index of each validation sample was 
then calculated. The measured and predicted melt indexes are shown in Table 8. 



WO 2005/049663 



43 



PCT/US2003/OI4565 



Acquisition of the validation spectra was interspersed with acquisition of the 
calibration spectra, at the corresponding "Vn M positions. 



Table 8: MI Validation, Reactor 2 



MI (Lab) 
(dg/min) 


MI (Model) 
(dg/TAin) 


MI (Lab) 
(dg/min) 


MI (Model) 
(dg/rnin) 


MI (Lab) 
(dg/min) 


MI (Model) 
(dg/min) 


VI: 0.733 


0.771 


17.291 


17.738 


23.390 


22.991 


0.754 


0.782 


17.896 


18.229 


V5: 1.907 


1.915 


0.798 


0.810 


20.620 


20.046 


1.908 


1.946 


0.727 


0.718 


V3: 52.180 


51.199 


1.958 


1.976 


0.721 


0.750 


52.020 


54.219 


1.902 


1.911 


V2: 17.649 


18.223 


V4: 24.880 


24.521 


1.930 


1.979 


18.399 


18.519 


20.760 


20.008 


1.930 


1.947 


19,844 


19.492 


18.667 


18.903 






21.480 


21.018 


16.682 


16.822 







5 - — — ■ 

[0153] Figure 8B depicts the data from Tables 7A, 7B and 8 graphically. The line 
in the Figure is the model prediction. The calculated R 2 value was 0.997, with a 
standard error of 2-86%. 



10 Examples 6-7 

[0154] Examples 6-7 demonstrate the effectiveness of the inventive methods on- 
line in a polymerization reaction system, for density determination. 
[01551 The measurements were carried out as described above in connection with 
Examples 4-5, except that a PCA/LWR model was developed for density. The 
15 samples used, and spectra acquired, are a subset of those of Examples 4-5. 
Laboratory measurements of density were made on the samples in addition to (he 
melt index measurements described above. 

Example 6: Density Model. Reactor 1 
20 [0156] One hundred forty-six polymer samples were evaluated The samples 
were divided into a group of 109 used for calibration (model development) and a 
group of 37 used for model validation. Each sample was a raetaUocene-catalyzed 
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LLDPE resin, in a density range of from about 0.912 to about 0.921 gfcra 3 . 
Raman spectra and laboratory density measurements were collected as described 
above. 

[0157] The lab values of density and the Raman spectra of the calibration data set 
5 were used to create a locally-weighted regression model for density, using 
principal component loadings and principal component scores. The measured 
densities and predicted densities are shown in Table 9. The deviations are not 
shown in the table, but are readily calculated from the tabulated data. The data are 
shown in the order taken (by column, within each table), to illustrate the 
10 effectiveness of the model under changing polymer conditions. A symbol "Vn" 
before an entry indicates that the nth set of validation spectra were taken before 
the marked entry, as shown by the corresponding notation in Table 10. 
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Table 9: Density (p, g/cm 3 ) Calibration, Reactor 1 



ofLatrt 

[ Y J K ± - ta - IJ / 


PfModoh 


P(Lab) 




riff ah\ 




A O^A^ 


A Qlti'i 

u.yzuj 


V-Jf IDC 


u.y i jy 


0.9151 


0.9150 


0 Q9A1 


A 01 QO 


A OKI 

w.7 ID 3 


A Ql KA 


A Ql CO 


o,y i jo 


0 01 fift 

l/a7 JLDD 


V.7JOO 


0 9161 


(I 0169 


A Ql 


A Ql 


0.9202 


0 9200 


0.9189 


0.9186 


^A 01 


a oi *;q 


0.9196 


0.9191 


0.9201 


0.9205 


n oi R6 

W.J' I ou 


A QIQfl 


0.9196 


0.9195 


0.9202 


0 9204 


0 Ql Rl 


A OI Rd 


0.9195 


0.9196 


0.9201 


0 9199 

v.7 J. 77 


ft 01 0^ 


A OI Q<\ 

u.y iyj 


0.9190 


0.9195 


0.9208 


0,9208 


A 0101 


A 01 04 

v. y i y*f 


0.9 192 


0.9187 


0.9201 




A 01 SI 




0.9195 






0 0166 

v.7 luv 


ft Ql 60 


A Ql 


0.9195 


0.9199 


0.9158 


0.9162 


A 01SO 


ft oi on 


0.9190 


0.9194 


0.9160 


0.9158 


Vi7JO X 


A OIRA 


0.9190 


0.9184 


0.9159 


0.9163 


0.9182 


0 01 R4 


0.9187 


0.9184 


0.9157 


0.9160 


0.9186 


0 01R5 


0-9187 


0,9190 


0.9157 


0.9156 


0.9188 


A 0105 


0.9188 


0.9190 


0.9159 


0.9157 


0.9186 


n oi si 


0.9196 


0.9192 


0.9157 


0 9156 


0 01 84 


A Olft* 


0.9196 


0,9195 


0.9161 


0.9156 


V4 fl 0104 


A 01 07 


0.9202 


0.9196 


0.9160 


0 9163 


0 91 RR 


A 91 S6 


0.9202 


0.9202 


0.9159 


0.9154 


0.9187 


A 01R5 


0.9207 


0.9204 


0,9155 


0.9158 


0.9180 


0 01 SA 


0.9200 


0.9199 


0.9149 


0.9146 


0 0144 

ITT 


A 0147 


0.9200 


0.9199 


0.9156 


0 9153 


0 017R 


0 011ft 


0.9197 


0 9200 


0.9164 


0.9162 


n oi to 


A 0117. 


0.9195 


0.9200 




0 0t 64 


n ons 


a om 


0.9195 


0.9193 


0-9162 


0 0164 

v.7l U*t 


n 0141 




V1 0.9195 


0.9199 


^0.9153 


0.9155 


0.9149 


0,9151 


0.9192 


0.9187 


0.9160 


0.9163 


0.9149 


0.9151 


0.9160 


0.9161 


0.9158 


0.9161 


0.9163 


0.9164 


0.9155 


0.9157 


0.9154 


0.9151 


0.9167 


0.9164 


0.9164 


0.9159 


0.9157 


0.9157 


0.9168 


0.9170 


0.9167 


0-9171 


0.9149 


0.9145 


0.9168 


0.9173 


0.9162 


0.9165 


0.9153 


0.9154 


0.9155 


0.9161 


0.9156 


0.9153 


0.9161 


0.9165 


0.9166 


0.9167 


0.9156 


0.9160 


0.9150 


0.9153 


0.9173 


0.9175 


0.9L62 


0.9165 


0.9155 


0.9152 






0.9159 


0,9162 


0.9150 


0.9151 
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[0158] The Raman spectra of the validation data $et were also collected, and new 
principal component scores were calculated from the validation spectra. Using the 
locally-weighted regression model, the density of each validation sample was then 
calculated. The measured and predicted densities are shown in Table 10. 
5 Acquisition of the validation spectra was interspersed with acquisition of the 
calibration spectra, at the corresponding "Vn" positions. 



Table 10: Density (p, stem 3 ) Validation, Reactor 1 



p(Lab) 


P(Model) 


p(Lab) 


p(Model) 


p(Lab) 


p(Model) 


VI; 0,9199 


0.9202 


0.9160 


0.9159 


0.9149 


0.9147 


0.9205 


0.9202 


0.9155 


0.9152 


V4: 0.9193 


0.9191 


0.9205 


0.9207 


0.9158 


0.9155 


0.9182 


0.9184 


0.9199 


0.9198 


0.9157 


0.9153 


0.9192 


0.9193 


0.9200 


0.9198 


0.9158 


0.9158 


0.9197 


0.9196 


0.9196 


0.9195 


0.9157 


0.9154 


0.9196 


0.9200 


0.9200 


0.9199 


0.9158 


0.9156 


0.9195 


0.9196 


0.9198 


0.9201 


0.9158 


0.9153 


0.9189 


0.9185 


0.9190 


0.9187 


V3: 0.9157 


0.9160 


0.9192 


0.9192 


0.9195 


0.9192 


0.9149 


0.9 ISO 


0.9198 


0.9197 


V2: 0,9159 


0.9161 


0.9168 


0.9168 


0.9192 


0.9193 


0.9188 


0.9188 


0.9149 


0.9146 






0.9159 


0.9161 


0.9150 


0.9153 







10 f0159] Figure 9A depicts the data from Tables 9 and 10 graphically. The line in 
the Figure is the model prediction. The calculated R 2 value was 0.978, with a 
standard error of 0.00028 g/cm 3 . 

Example 7: Density Model, Reactor 2 
15 [0160] The procedure described in Example 6 was followed, except as noted, 
sampling this time from the Reactor 2 polymer. One hundred sixty-four polymer 
samples were evaluated. The samples were divided into a group of 153 used for 
calibration (model development) and a group of 13 used for model validation. 
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Each sample was a Ziegler-Natta-catalyzed LLDPE resin, in a density range of 
from about 0.916 to about 0.927 g/cm 3 . Raman spectra and laboratory density 
measurements were collected as described above. 

[0161] The lab values of density and the Raman spectra of the calibration data set 
5 were used to create a locally-weighted regression model for density, using 
principal component loadings and principal component scores. The measured 
densities and predicted densities are shown in Tables 11A-11B. The deviations 
are not shown in the table, but are readily calculated from the tabulated data. The 
data axe shown in the order taken (by column, within each table), to illustrate the 
10 effectiveness of the model under changing polymer conditions. A symbol "Vn" 
before an entry indicates that the nth set of validation spectra were taken before 
the marked entry, as shown by the corresponding notation in Table 12. Table 1 IB 
is a continuation of Table 11 A. 
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Table 1 1 A: Density (p, g/W) Calibration, Reactor 2 



p(Lab) 


P(Model) 


ofLsib) 


oflVloden 


p\J~>aX>) 


p(Moael) 


0.9182 


0.9182 


V2 0-9269 


/ KM 




0.9184 


0.9180 


0.9184 


0.9267 


0 9268 




0,9178 


0.9207 


0.9209 


0.9259 


0.9263 




Am id 

U.9179 


0.9220 


0.9225 


0.9246 


0.9249 




n ni on 


0.9220 


0.9221 


0.9235 


0.9235 


0 917R 

\S<-7 I/O 




0.9220 


0.9217 


0.9246 


0.9250 


0 9178 


A Q1 TO 


1 0.9218 


0.9219 


0.9248 


0.9246 


0 9190 




0.9218 


0.9219 


0.9256 


0.9257 


0.9197 


ft 01 07 
kj*yly<& j 


0.9217 


0.9217 


0.9251 


0.9253 


0.9184 


ft Q17K 


0.9220 


0.9220 


0.9246 


0.9246 


0.9184 


0 QIQft 


0.9226 


0.9221 


0.9253 


0.9253 


0.9199 


ft 91 QR 

Vi7170 


0.9217 


0.9216 


0,9260 


0.9259 


0.9182 


0 9186 

Ui7I OkJ 


0.9219 


0.9222 


0.9265 


0.9268 


0.9177 


ft 0174 


0,9225 


0.9224 


0.9265 


0.9264 


0 9180 

V/ ■ J luv 


A Q17R 


0.9221 


0.9223 


0.9261 


0.9260 


v "0 9159 


0 91 £1 


0.9216 


0.9217 


0.9251 


0.9256 


0 9169 


\J.2r 1U7 


0.9218 


0.9219 


0.9252 


0.9254 


0 9173 


ft 01 77 


0.9218 


0.9213 


0.9249 


0.9254 


0 9172 


ft Q17fi 


0.9216 


0.9220 


0.9257 


0.9255 


0 9178 


v«y lot 


Vl 0,9268 


0.9262 


0.9252 


0.9247 


0,9181 




0.9256 


0.9259 


0.9244 


0.9249 


0.9179 


ft Q17d 


0.9254 


0.9254 


0.9245 


0.9245 


0.9204 


0 Q704. 


0.9255 


0.9257 


0.9250 


0.9248 


0.9174 


0,9180 


0.9252 


0.9253 


0.9256 


0.9260 


0.9183 


0.9185 


0.9247 


0.9253 


0.9256 


0.9260 


0.9184 


0.9180 


0.9255 


0.9260 


^0.9216 


0.9213 


0,9177 


0.9176 


0.9250 


0.9246 


0.9168 


0.9153 


0.9173 


0.9169 


0.9264 


0.9267 


0.9184 


0.9186 


0.9176 


0.9174 


0.9259 


0.9258 


0.9191 


0.9189 


0.9178 


0.9181 


0.9253 


0.9250 


0.9188 


0,9187 


0.9180 


0.9181 


0.9247 


0.9245 


0.9185 


0,9186 


0.9182 


0.9185 
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Table 1 IB: Density (p, g/cm 3 ) Calibration, Reactor 2, continued 



p(Lab) 


P(Model) 


pCLab) 


p(Modei) 


p(Lab) 


p(Model) 


0.9182 


0.9185 


0.9207 


0.9203 


0.9253 


0.9253 


0.9166 


0.9170 


0.9213 


0.9213 


0.9253 


0.9253 


0.9187 


0.9185 


0.9218 


0,9215 


0.9265 


0.9269 


0.9184 


0.9189 


0.9226 


0.9231 


0,9261 


0.9263 


0.9181 


0.9183 


0.9221 


0,9225 


0.9259 


0.9257 


0.9182 


0.9177 


0.9218 


0.9217 


0.9257 


0.9257 


0.9181 


0.9180 


0.9216 


0.9217 


0.9260 


0.9255 


0,9185 


0.9182 


0.9223 


0.9218 


0.9251 


0.9252 


0,9180 


0.9183 


0.9223 


0.9220 


0.9238 


0,9236 


0,9182 


0.9187 


0.9222 


0.9227 


0.9248 


0.9243 


0.9183 


0.9187 


0,9220 


0.9218 


0.9270 


0.9275 


0.9183 


0.9188 


0.9223 


0.9224 


0.9247 


0.9243 


0.9180 


0.9185 


0.9222 


0.9223 


0.9244 


0.9241 


0.9181 


0.9177 


0.9220 


0.9218 


0.9244 


0.9239 


0.9180 


0.9183 


0.9224 


0.9223 


0.9249 


0.9245 


0.9182 


0.9182 


0.9224 


0.9221 


0.9250 


0.9246 


0.9179 


0.9177 


0.9225 


0,9223 


0.9248 


0.9246 


0.9182 


0.9178 


0.9223 


0.9220 


0.9250 


0.9256 


0.9181 


0.9182 


0.9216 


0.9216 






0.9204 


0.9202 


0.9253 


0.9253 







[0162] The Raman spectra of the validation data set were also collected, and new 
5 principal component scores were calculated from the validation spectra. Using the 
locally-weighted regression model, the melt index of each validation sample was 
then calculated. The measured and predicted melt indexes are shown in Table 12- 
Acquisition of the validation spectra was interspersed with acquisition of the 
calibration spectra, at the corresponding " Vn" positions. 
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Table 12: D&nshy (p, g/cm 3 ) Validation, Reactor 2 



p(Lab) 


P(Model) 


p(Lab) 


p(ModaJ) 


p(Lab) 


p(ModeI) 


VI; 0,9216 


0.9217 


0.9262 


0.9265 


0.9184 


0.9181 


0.9221 


0.9219 


0.9250 


0.9256 


0.9185 


0.9183 


0.9220 


0.9218 


V3: 0.9238 


0.9243 


0.9185 


0.9188 


V2: 0.9254 


0.9257 


0.9228 


0.9230 






0.9250 


0.9247 


V4: 0.9182 


0.9180 







[0163J Figure 9B depicts the data from Tables 11 A, I IB and 12 graphically. The 
line in the Figure is the model prediction, The calculated R 2 value was 0.989, 
5 with a standard error of 0.00034 g/cm 3 . 
Examples 8-9 

[0164] Examples 8-9 demonstrate the effectiveness, precision and accuracy of 
processes of the invention to predict melt index and density on-line, in a 
commercial-scale fluidised-bed polymerization reactor. The Raman system was 
1 0 as described above but used a 400 mW diode laser operating at 785 nm. The fiber 
optic cable used to couple the electrical components of the iiistrument to the 
Raman probe (approximately 150 m distant) was a 62 nm excitation/1 OO^m 
collection step index silica fiber. 

[0165] Melt index and density models were developed by continuously collecting, 
15 and saving Raman data as individual spectra every 3-10 minutes, on each of two 
reactors. Validation of each model was accomplished by then using the model on- 
line to determine the polymer properties. 
Example 8 

[0166] Polymer melt index was predicted on-line in a commercial-scale fluidized- 
20 bed reactor forming various grades of polyethylene copolymer. The prediction 
was carried out approximately every 12 minutes for about 5 weeks, Nearly 500 
samples were also tested the laboratory, using the standard ASTM D-1238, 
condition E (2.16 kg load, 190 °C) protocol. The results, are shown in Table 13, 
where "MI model" indicates the melt index 12.16 predicted by the model* and "MI 
25 lab" indicates the value obtained in the laboratory by the ASTM method. The 
same data are shown graphically in Figure 10, except that the Figure also shows 
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Hie predicted MI for samples not corresponding to lab measurements. The 
predicted MI values are spaced sufficiently closely in time that they appear in the 
Figure to bo a line. 

Table 13 

Jj™ **} MIbb S MI MI lab Time MI MI lab time m mu^T 

{ ^ iSStx ***** (dflyS) Wlttin) (days > < d ^) (*55 model (SgK) 
aSBSg (dfi /mm) (<lfefoHn) (dg/mjp) 

0.009 0.958 0,957 16.836 1.099 1330 ESS U03 ToSS 27.837 2 10& 2 090 

awe low 1.020 laoas i.uj i.uo 19.126 1,125 i.oso 27.90s 2 we zoeo 

0.155 1.036 1,030 10.999 1.037 1.000 19.169 1.077 UOO 28.009 2153 2350 

0.258 1.007 0.998 11.042 1.107 L085 19.235 1.030 1.060 28.086 2 136 2 140 

0.327 1.002 0.996 11.085 r.314 1.302 19.341 1.019 1,010 28 J 64 2 096 M35 

6.404 0.982 1.006 11.120 1.451 3.426 19393 0.997 1.050 28.241 2 125 2340 

0.4^9 0.941 0.927 11.163 1.497 1.495 1930S 0.926 0.927 28,293 1 861 3360 

0.585 0.944 0.940 31.249 1.444 1.456 19339 0,875 0.886 28.336 3 483 1490 

1.620 



0.654 1.017 1.01 J 11.326 1.374 1.255 19.582 0.605 0.604 28.362 U73 
0.748 1.033 J.037 JI395 M*2 1.468 19.617 0.563 0342 28.405 1.472 1470 

2™ JUS ™ll I!*!?! \ A A* U40 l9 - 660 0,432 0Jto 1.482 1300 

1.470 
1 .320 
1.370 
1.440 



0.903 0.940 0.970 11.5*4 1383 MOO 19.746 0.550 0.551 28 508 I486 

0.998 0553 0.930 U.670 1.364 1J70 19.832 0375 0.581 28.542 1332 

1.084 0.944 04SO 11.748 1.332 1330 19.900 0362 0345 28485 1375 

1-161 0.971 0.980 11.834 1.370 1.370 19.995 0341 0J39 28*71 1427 1W 

!!"!?? i' 6 ! 5 I,6fi0 20i09D 0374 28 - w 1404 mio 

* ~ 1310 

1.380 



1.247 1^20 1.210 113U 2443 2.43 1 20.167 0493 0.572 28.835 3313 

1.290 1.373 1.360 11.963 3.017 3.029 20.233 0.6ZO 0385 28.903 1380 w 

1.324 1.461 L460 11.997 3.102 3.107 20.330 0.526 0330 38.998 1331 1340 

1.367 1.502 1.520 12.040 3307 3 201 20.399 0300 0.506 29.084 1367 U60 

1.419 1.486 1306 12.083 3.357 3.350 20.5 U 0357 0,529 29-170 1307 1320 

1.995 2.763 2.770 12.126 3.323 3348 20.580 0.506 0.511 29.247 1 319 1320 

2.081 2,650 2.650 12.169 3330 3.333 20.657 0.538 0.526 29333 1.374 3 370 

2.159 2.766 2.810 12.246 3305 3.482 20.752 0333 0.526 29-394 1 261 1250 

2.236 2.777 2.590 12332 3491 3.211 20.838 0316 0316 29305 1.216 1220 

2339 2.738 2.740 12.401 3.691 3.691 20.898 0331 0336 29374 1229 1 220 

2.399 2.799 2.800 12304 3.713 3.660 21-000 0 309 0.543 29.669 1,227 1 230 

2311 2.970 2.971 12.582 4.040 4.080 21.095 0.523 0.522 29 755 1.215 1230 

2346 3.056 3.072 12.668 3.947 3.960 21.163 0394 0.542 29.832 1 206 1210 

2380 3.233 3.226 12.7S4 3.775 3.770 21.249 0.479 0.530 29.901 1322 1*220 

2*66 3.427 3.423 12.805 3.743 3.730 21327 0315 0320 29.987 1.205 1 210 

2.752 3 342 3.545 !Z$40 3.698 3.600 21.404 0.960 0.920 30.082 1 347 1350 

2.838 3.619 3.699 12.900 3.753 3,770 21-456 L040 1.056 30.125 1.221 1340 

2.907 3393 3.5B0 12.960 3.964 3-950 21499 1,150 1.133 30,168 1 177 1 180 

3.001 3.446 3380 13,003 3.693 3.150 21342 1.152 1.172 30.211 1.0* 1080 

3.079 3.5)4 3 300 13.029 3 301 3310 21385 1.157 1.080 30.254 1 047 lOlO 

3.165 3.70S 3.710 13.081 3.112 3.110 21.671 1.005 1.124 30.322 1057 1070 

3.242 3.702 3.710 13.132 3,063 2.870 21.748 1.030 1.067 30.400 1 012 1 020 

3.337 3.710 3.710 13.167 3.523 3.520 21.834 1.060 1.036 30303 0-987 1 012 

3.397 3.713 3.710 13.210 3.624 3.630 21.894 1.090 1.093 30380 0.943 0.932 

3309 3.478 3.482 13.253 3 380 3.600 21.998 1.060 1.1 1Q 30.666 0,889 0902 

\™ 5361 13.296 3.705 3.720 22.041 1309 1,230 31.707 3.260 1348 

2.763 



3364 3466 3.465 13330 3.700 3.700 22.084 1,548 1.650 31.750 2J76 

3.750 3,430 3.443 13 459 3.222 3.220 22.135 2.062 2410 31784 3394 

3.836 3458 3.4S9 13302 3.1B8 3.180 22.170 2.106 2.U0 31.836 3350 

3.913 3306 3300 13343 3334 3-240 22.221 2.066 2.120 31.879 4310 

3.999 3301 3390 13379 3.256 3.250 22.247 2.108 2.050 31.903 3 968 

4.042 3,411 3320 13.622 3313 3390 22.324 2.117 2,080 31.956 4140 4 110 
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4.085 3.466 3.460 13.674 3.327 3.350 22.402 2.155 2.320 32.017 4.019 4.040 

4.171 3.751 3.730 13.751 3.404 3350 22.496 2.156 2.190 32.085 4.431 4.440 

4.248 3.713 3.700 13.837 3.488 3.550 22.582 2.013 2.040 32.171 4^645 4*650 

4.334 3.493 3.500 13.898 3.421 3.390 22.660 . 2.067 2.080 32.249 4.737 4.730 

4.403 3.497 3.460 14.001 3.425 3.440 22.754 2.108 2.100 32.326 4*827 4.840 

4.498 3.459 3.406 14.087 3.487 3.500 22.823 2.108 2.120 32.404 4.756 4.754 

4.584 3.516 3.528 14.164 3.459 3.480 22.858 2.348 2350 32.498 4.212 4*182 

4.670 3.544 3.555 14.250 3.429 3.420 22.901 3.418 3.380 32.576 3.953 3.975 

4.747 3.601 3.616 14.345 3.373 3.380 22.961 3.794 4.122 32.627 4.189 4^217 

4.833 3.612 3.593 14.396 3.363 3.370 23.004 3.507 3.588 32.670 4.309 4^297 

4.902 3.514 3.536 14.500 3.420 3.259 23.038 3.358 3316 32.748 4.328 4.320 

4.997 3.560 3.559 14.586 3.467 3.462 23.081 3.234 3.192 32.825 4.300 4315 

5.091 3.584 3.677 14.663 3.566 3.400 23.133 3.187 3.196 32.902 4.338 4.366 

5.177 3.583 3.350 14.749 3.478 3.475 23.167 3.424 3.417 32.997 4.263 4.270 

5.255 3.531 3.554 14.826 3.383 3.409 23.253 3.507 3.497 33.049 4.225 4.230 

5.332 3.473 3.476 14.904 3.323 3.341 23.331 3.587 3.581 33.083 3.947 3^930 

5.409 3.528 3.521 14.998 3.636 3.640 23.408 3.353 3.351 33.126 3.605 3.610 

5.504 3.537 3.526 15.084 3.537 3.550 23.503 3.292 3.300 33.169 3.446 3.460 

5.581 3.484 3.436 15.170 3.369 3.360 23.580 3.344 3.330 33.212 3.416 3.400 

5.667 3.532 3.592 15.256 3.270 3.300 23.666 3.395 3.380 33.246 3.511 3.530 

5.753 3.482 3.516 15.334 3.644 3.630 23.752 3.396 3.390 33.332 3.503 3.520 

5.839 3.505 3.516 15.394 3.289 3.310 23.838 3.406 3.270 33.401 3.387 3.380 

5.900 3.516 3.464 15.497 3.136 3.150 23.898 3.386 3.390 33.496 3.327 3.340 

6.003 3.455 3.464 15.583 3.439 3.435 23.993 3.509 3.478 33.590 3303 3.292 

6.080 3.466 3.394 15.661 3.460 3.448 24.079 3.576 3.586 33.668 3.472 3.457 

6.166 3.468 3.439 15.755 3.461 3.494 24.165 3.530 3.401 33.754 3.625 3.595 

6.252 3.666 3.655 15.833 3.466 3.447 24.242 3.513 3.520 33.831 3.544 3.520 

6.330 3.723 3.714 15.901 3.612 3.620 24.320 3.417 3.414 33.900 3.603 3.626 

6.407 3.712 3.778 16.005 3.458 3.450 24.397 3.463 3.439 33.986 3.551 3.570 

6.510 3.578 3.555 16.039 3.301 3.310 24.500 3.406 3.410 34.081 3.618 3.610 

6.588 3.483 3.480 16.082 3.222 3.220 24.586 3.399 3.400 34.158 3.478 3.470 

6.665 3.454 3.470 16.125 2.997 2.980 24.664 3.241 3.290 34.253 3,572 3.580 

6.751 3.399 3.540 16.168 2.804 2.790 24.750 3.272 3.280 34.321 3.576 3.560 

6.854 3.379 3.360 16.211 2.751 2.740 24.836 3.322 3.280 34.399 3.663 3.660 

6,906 1.824 1.810 16.254 2372 2.390 24.896 3.480 3.481 34.493 3.596 3.609 

6.949 1.157 1.140 16.288 2360 2.360 24.999 3.345 3.313 34.588 3.312 3331 

7.000 0.899 0.860 16331 2387 2.400 25.076 3.316 3,308 34.665 3.238 3.261 

7.043 0.888 0.883 16,400 2.467 2350 25.162 3.479 3.479 34.751 3.358 3362 

7,086 0.937 0.940 16.495 2.569 2.564 25.240 3.491 3.509 34.820 3.424 3.416 

7.129 1.001 0.980 16,581 2.609 2.610 25.343 3.544 3.525 34.906 3.465 3.458 

7.164 1.053 1.060 16.667 2.680 2.660 25395 3.373 3.405 35.009 3.444 3.440 

7.250 1.105 1.110 16.710 2.649 2.638 25.498 3.582 3.580 35.087 3.450 3.470 

7336 1.083 1.070 16.744 2.074 2.079 25.584 3 371 3370 35.164 3.474 3.470 

7.413 1.051 1.040^ 16.796 1.916 1.938 25.670 3.223 3.190 35.250 3.485 3.510 

7.508 1.003 1.010 16.830 1.917 1.995 25.747 3.274 3 300 35.328 3.627 3.630 

7.585 0.975 0.973 16.865 1.960 1.967 25.833 3.518 3.500 35396 3.618 3.620 

7.680 0.958 0.950 16.899 2.049 2.070 25.919 3.414 3.409 35.500 3,654 3.669 

7.749 0.999 0.960 16.994 2.211 2.200 25.997 3.445 3.451 35.586 3.354 3311 

7.835 1.026 1.000 17.080 2.152 2.150 26.074 3.460 3.488 35.663 3.389 3.404 

7.903 1.004 1.000 17.166 2322 2 320 26.169 3.591 3.603 35.749 3.463 3.452 

7.998 1.022 1.033 17.252 2 3 52 2.240 26.246 3.825 3.836 35.835 3.550 3.544 . 

8.101 1.038 1.037 17329 2 342 2340 26.332 3.681 3.668 35.895 3.449 3.484 

8.161 1.021 1.025 17.398 2.239 2.247 26.401 3.392 3372 36.007 3371 3380 

8.265 0.968 0.990 17.501 2.110 2.112 26.487 3.395 3380 36.084 3.382 3.390 

8325 0.992 0.988 17,587 2.081 2.080 26.573 3.404 3.390 36.170 3.448 3.440 

9.460 0.812 0.754 17.664 2.098 2.119 26.667 3.320 3.300 36.248 3.634 3.630 

9.503 0.787 0.790 17.750 2.147 2.130 26.753 3370 3370 36334 3.743 3.730 

9.546 0.886 0.871 17.836 2.273 2301 26.831 3.428 3.410 36.403 3.633 3.630 
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10 



15 



20 



9.580 
9.666 
9.752 
9.830 
9.899 
10.002 
10.088 
10.165 
10.243 
10.329 
10.397 
10.501 
10.578 
10.664 
10.750 



0.992 

0.996 

1.012 

1.006 

1.024 

1.127 

1.016 

1.051 

1.062 

1.130 

1.052 

1.057 

1.131 

1.130 

1.141 





17.905 


2.202 


2.205 


26.900 


3.539 


3.280 


36.506 




18.000 


2.058 


2.050 


27.003 


3.731 


3.723 


36.583 


1.020 


18.086 


2.002 


2.080 


27.080 


3.498 


3.476 


36.669 


1.020 


18.163 


2.043 


2.040 


27.158 


3.572 


3.552 


36.747 


1.010 


18.249 


2.094 


2.110 


27.252 


3.349 


3.361 


36.781 


1.110 


18.335 


2.203 


2.180 


27.287 


2.622 


3.265 


36.815 


0.990 


18.395 


2.260 


2.250 


27.330 


2.483 


2.481 


36.858 


1.030 


18.490 


2.264 


2.257 


27.381 


2.497 


2.505 


36.919 


1.080 


18.662 


2.292 


2.277 


27.407 


2.433 


2.438 


36.962 


1.130 


18.757 


2.195 


2.158 


27.467 


2.140 


2.170 


37.048 


1.010 


18.843 


2.208 


2.189 


27.510 


2.047 


2.050 


37.125 


1.040 


18.903 


2.112 


2.139 


27.536 


1.920 


1.920 


37.220 


1.130 


18.963 


1.676 


1.690 


27.579 


1.883 


1.880 


37.288 


1.140 


18.997 


1.390 


1.400- 


27.665 


2.051 


2.030 


37.366 


1.100 


19.040 


1.219 


1.240 


27.751 


2.052 


2.060 



3.376 

3.399 

3.300 

1.486 

1.435 

1.283 

1.293 

1.394 

1.404 

1.417 

1.457 

1.551 

1.571 

1.571 



3.382 
3398 
3314 
1.483 
1.429 
1.296 
1.306 
1390 
1.420 
1.430 
1.480 
1.570 
1.570 
1.540 



[0167] Table 13 and Figure 10 show the accuracy and precision of the on-line 
process over a long period of time, and a range of melt index values. The gaps in 
the Figure indicate periods when the reactor was down. The horizontal regions 
indicate continued production of a particular grade, and the steep vertical regions 
correspond to transitions between different grades. The data further show that the 
inventive on-line processes are accurate and precise even during grade transitions. 
The 3a accuracy of the predictions relative to the lab values over the entire 5- 
week period was ± 0.069 g/1 0 min, 

[0168] Additionally, to test for model precision and long-term drift, the predicted 
MI of approximately 2200 samples of a particular grade was monitored for a static 
sample over a four-week period, in each of two commercial-scale fluidized bed 
reactors. In each reactor, the data showed a 3o standard deviation of 0.012 g/10 
min (for sample with melt indexes of 1.0 and 0.98 g/10 min; i.e., about 1%), and 
no measurable long-term drift. 
Example 9 

[0169] Polymer density was predicted on-line along with the melt index 
predictions of Example 8, applying a density model to the same samples and 
spectra as in Example 8. Nearly 300 samples were also tested the laboratory, 
using the standard ASTM D1505 and ASTM D1928, procedure C protocol. The 
results, are shown in Table 14, where "p model" indicates the density predicted by 
the model, and "p lab" indicates the value obtained in the laboratory by the ASTM 
method. The same data are shown graphically in Figure 1 1, except that the Figure 
also shows the predicted density for samples not corresponding to lab 
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measurements. The predicted density values are spaced sufficiently closely in 
time that they appear in the Figure to be a line. 

Table 14 



Time p model p lab Time p model plab Time pmode] plab Time p model p lab 
(days) (g/cm 3 ) (g/cm 3 ) (days) (g/cm 3 ) (g/cm 3 ) (days) (g/cm 3 ) (g/cm 3 ) (days) (g/cm 3 ) (g/cm 3 ) 



0.009 


0.9173 


0.9173 


10.664 


0.9170 


0.917 


19.505 


0.9168 


05168 


28.164 


0.9168 


05169 


0.155 


0.9183 


0.9182 


10.836 


0.9171 


0.9171 


19.582 


0.9199 


05199 


28.241 


05174 


0.9173 


0.327 


05176 


0.9175 


10.999 


0.9186 


0.9186 


19.660 


0.9214 


0.9214 


28.336 


0.9167 


0.9168 


0.499 


0.9175 


0.9175 


11.085 


0.9175 


0.9175 


19.746 


0.9202 


0.9204 


28.405 


05161 


05160 


0.654 


0.9172 


0.9173 


11.163 


0.9177 


0.9177 


19.832 


0.9199 


05200 


28.508 


05164 


05164 


0.826 


0.9178 


0.9178 


11.249 


0.9179 


0.9179 


19.995 


0.9206 


0.9206 


28.671 


0.9169 


05168 


0.998 


0.9173 


05173 


11326 


0.9184 


0.9183 


20.167 


0.9208 


05207 


28,835 


05168 


05168 


1.161 


0.9173 


0.9173 


11.507 


0.9176 


0.9177 


20330 


0.9208 


0.9206 


28.998 


0.9164 


05164 


1.247 


0.9167 


0.9166 


11.670 


0.9175 


0.9173 


20.511 


0.9207 


0,9207 


29.170 


0.9167 


05163 


1.324 


0.9169 


0.9169 


11.834 


0.9173 


0.9172 


20.657 


0.9203 


0.9203 


29333 


05169 


0.9170 


1.995 


0.9172 


0.9172 


11.911 


0.9176 


0.9175 


20.838 


0.9215 


05214 


29.505 


05164 


0.9163 


2.159 


0.9168 


0.9168 


11.997 


0.9173 


0.9173 


21.000 


0.9206 


05205 


29.669 


05171 


05170 


2.339 


0.9169 


0.9168 


12.083 


0.9180 


0.9182 


21.163 


0.9208 


0.9207 


29.832 


05173 


05 1 73 


2,511 


0.9187 


0.9186 


12.169 


0.9181 


0.9182 


21327 


0.9211 


0.9210 


29587 


0.91 74 


05177 


2.580 


0.9185 


0.9184 


12332 


0.9186 


0.9185 


21.404 


0.9186 


0.9188 


30,168 


05165 


05164 


2.666 


0.9183 


0.9184 


12.504 


0.9172 


0.9172 


21.499 


0.9170 


0.9168 


30.254 


0.9172 


05172 


2.838 


0.9179 


0.9180 


12.668 


0.9167 


0.9166 


21.585 


0.9168 


0.9167 


30322 


0.9170 


05171 


3.001 


0.9166 


0.9167 


12.840 


0.9165 


0,9166 


21.671 


0.9172 


0.9172 


30.400 


05162 


05162 


3.165 


0.9173 


0.9173 


13.003 


0.9173 


0.9173 


21.834 


0.9170 


0.9170 


30.503 


05171 


05173 


3337 


0.9172 


0.9172 


13.167 


0.9176 


0.9176 


21.998 


0.9171 


0.9171 


30.666 


0.9181 


05180 


3.509 


0.9181 


0.9181 


13330 


0.9176 


0.9175 


22.084 


0.9174 


' 0.9174 


31.750 


05205 


05205 


3.664 


0.9173 


0.9173 


13.502 


0.9174 


0.9172 


22.170 


05164 


05164 


31.836 


05195 


05195 


3.836 


0.9173 


0.9172 


13.674 


0.9173 


0.9174 


22.247 


05167 


0.9167 


31.905 


0.9189 


05188 


3.999 


0.9165 


0.9165 


13.837 


0.9176 


0.9176 


22324 


05168 


05167 


32.017 


0.9174 


05176 


4.171 


0.9176 


0.9177 


14.001 


0.9176 


0.9175 


22.496 


05176 


0.9177 


32.171 


05176 


05177 


4334 


0.9175 


0.9173 


14.164 


0.9174 


0.9175 


22.660 


0.9166 


05167 


32326 


05176 


0.9175 


4.498 


0.9171 


0.9172 


14.345 


0.9172 


0.9170 


22.823 


05168 


0.9168 


32.498 


05161 


0.9160 


4.670 


0.9177 


0.9175 


14.500 


0,9173 


0.9173 


22.901 


05169 


05168 


32.670 


05171 


05171 


4.833 


0.9179 


0.9179 


14.663 


0.9178 


0.9179 


23.004 


05175 


0.9175 


32.825 


0.9175 


05174 


4.997 


0.9179 


0.9178 


14.826 


0.9185 


0.9183 


23.167 


05184 


05184 


32.997 


05171 


05171 


5.177 


0.9175 


0.9176 


14.998 


0.9174 


0.9173 


23331 


0.9180 


0.9178 


33.083 


0.9170 


05169 


5.332 


0.9172 


0.9173 


15.170 


0.9172 


0.9171 


23.503 


05177 


05178 


33.169 


0.9171 


0.9170 


5J04 


0.9177 


0.9177 


15.334 


0.9171 


0.9171 


23.666 


05175 


0.9175 


33.246 


05170 


05170 


5.667 


0.9173 


0.9173 


15.497 


0.9174 


0.9173 


23.838 


05174 


0.9175 


33332 


0.9170 


0.9170 


5.839 


0.9171 


0.9171 


15.661 


0.9170 


0.9171 


23.993 


0.9182 


0.9183 


33.496 


0.9175 


0.9175 


6.003 


0.9166 


0.9166 


15.833 


0.9171 


0.9171 


24.1.65 


05184 


05184 


33.668 


05174 


0.9176 


6.166 


0.9169 


0.9169 


16.005 


0.9174 


0.9175 


24320 


05172 


0.9172 


33.831 


05168 


05170 


6.330 


0.9179 


0.9180 


16.082 


0.9171 


0.9172 


24.500 


05175 


0.9173 


33.986 


0.9169 


05168 


6.510 


0.9175 


0.9175 


16.168 


0.9176 


0.9175 


24.664 


05178 


0.9178 


34.158 


05171 


0.9170 


6.665 


0.9177 


0.9177 


16.254 


0.9181 


0,9181 


24.836 


05185 


0.9184 


34321 


05174 


0.9175 


6.854 


0.9171 


0.9170 


16331 


0.9180 


0.9179 


24.999 


05180 


0.9181 


34.493 


05169 


05170 


6.949 


0.9157 


0.9157 


16.495 


0.9171 


0.9171 


25.162 


0.9172 


0.9172 


34.665 


05170 


0.9170 


7.000 


0.9165 


0.9165 


16.667 


0.9171 


0.9169 


25.343 


0.9176 


0.9175 


34.820 


05172 


0.9171 


7.086 


0.9167 


0.9167 


16.744 


0.9171 


0.9169 


25.498 


0.9169 


0.9170 


35.009 


05172 


05173 


7.164 


0.9174 


0.9374 


16.830 


0.9163 


0.9163 


25.670 


05173 


0.9173 


35.164 


0.9176 


0.9176 


7.336 


0.9179 


0.9180 


16594 


0.9164 


0.9164 


25.833 


05171 


05171 


35328 


05177 


0.9176 


7.508 


0.9182 


0.9181 


17.166 


0.9165 


0.9163 


25.997 


05172 


0.9172 


35.500 


05176 


0.9176 


7.680 


0.9179 


0.9178 


17329 


0.9162 


0.9161 


26.169 


0.9173 


0.9173 


35.663 


05183 


0.9182 


7.835 


0.9178 


0.9178 


17.501 


0.9164 


0.9164 


26332 


05172 


0.9172 


35.835 


0.9169 


0.9168 


7.998 


0.9175 


0.9176 


17.664 


0.9169 


0.9169 


26.487 


0.9172 


0.9173 


36.007 


0.9166 


05164 


8.161 


0.9173 


0.9174 


17.836 


0.9165 


0.9367 


26.667 


05173 


0.9174 


36.170 


05174 


05174 
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8 - 325 0.9168 0.9169 18.000 0,9175 0.9173 26.831 0.9173 0.9)73~ 36J34 0 9169 09171 

9.503 0.9188 0.9190 18.163 0.9168 0.9168 27.003 0.9174 0.9174 36.506 0 9172 0 9173 

9J80 0.9185 0.9185 18.335 0.9170 0.9171 27.158 0.9164 0.9164 36.669 0 9169 0 9169 

9.666 0.9179 0.9181 18.490 0.9168 0.9168 27.330 0.9173 0.9174 36.747 0 9162 0 9162 

9.752 0.9175 0.9174 18.662 0.9176 0.9176 27.407 0.9162 0.9162 36 815 0 9162 0 9162 

9.830 0.9175 0.9174 18.843 0.9172 0.9171 27.510 0.9162 0.9160 36 962 0^9172 0 9172 

10.002 0.9173 0.9174 18.997 0.9181 0.9181 27.579 0.9169 0.9169 37.125 0 9174 0 9172 

10.165 0.9173 0.9171 19.083 0.9176 0.9174 27.665 0.9169 0.9168 37 220 0 9175 0 9174 

10.329 0.9172 0.9173 19.169 0.9163 0.9164 27.837 0.9169 0.9168 

10.501 0.9172 1 0.9173 19.341 0.9163 0.9163 28.009 0.9169 0.9171 

[0170] Table 14 and Figure 1 1 show the accuracy and precision of the on-line 
process over a long period of time, and a range of density values. As in the 
previous Example, the gaps in the Figure indicate periods when the reactor was 
5 down, the horizontal regions indicate continued production of a particular grade, 
and the steep vertical regions correspond to transitions between different grades. 
The data further show that the inventive on-line processes are accurate and precise 
even during grade transitions. The 3a accuracy of the predictions relative to the 
lab values over the entire 5-week period was ± 0.00063 g/cm 3 . 
10 [0171] Additionally, to test for model precision and long-term drift, the predicted 
density of the same approximately 2200 samples of Example 8 was monitored for 
a static sample over a four-week period, in each of two commercial-scale fluidized 
bed reactors. In each reactor, the data showed a 3a standard deviation of 0.00006 
g/cm 3 (for samples with densities of 0.9177 and 0.9178 g/cm 3 ), and no measurable 
15 long-term drift. 

[0172] Having thus described the invention with reference to specific examples, 
the following is intended to set forth particular preferred embodiments, without 
intending to limit the spirit and scope of the appended claims. Although described 
below with reference to in situ sampling, the descriptions below also apply to the 
20 extractive sampling except where it would be readily apparent to one of ordinary 
skill in the art in possession of the present disclosure that extractive sampling 
would not apply. 

[0173] One preferred embodiment is. a process for deterrnimng polymer properties 
in a polymerization reactor system, the process comprising: (a) obtaining a 
25 regression model for determining a polymer property, the regression model 
including principal component loadings and principal component scores; (b) 
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acquiring a Raman spectrum of a sample comprising polyolefin; (c) calculating a 
new principal component score from at least a portion of the Raman spectrum and 
the principal component loadings; and (d) calculating the polymer property by 
applying the new principal component score to the regression model. Even more 
5 preferred emodiments include one or more of the following: wherein the step of 
obtaining a regression model comprises: (i) obtaining a plurality of Raman spectra 
of samples comprising polyolefms; (ii) calculating principal component loadings 
and principal component scores from the spectra obtained in (i) using principal 
component analysis (PCA); and (iii) forming the regression model using the 

10 principal component scores calculated in (ii) such that the regression model 
correlates the polymer property to the principal component scores; wherein the 
regression model is a locally weighted regression model; wherein the polymer 
property is selected from density, melt flow rate, molecular weight, molecular 
weight distribution, and functions thereof; wherein the sample comprises 

15 polyolefin particles; wherein the step of acquiring a Raman spectrum comprises: 
(i) providing the sample of polyolefin particles; and (ii) irradiating the sample and 
collecting scattered radiation during a sampling interval using a sampling probe, 
wherein there is relative motion between the sample and the sampling probe 
during at least a portion of the sampling interval; wherein the polymerization 

20 reactor is a fluidized-bed reactor; wherein the reactor includes a cyclone; wherein 
the process further comprises: (i) obtaining a second regression model for 
determining a second polymer property, the second regression model including 
second principal component loadings and second principal component scores; (ii) 
calculating a new second principal component score from at least a portion of the 

25 Raman spectrum and the second principal component loadings; and (iii) 
calculating the second polymer property by applying the new second principal 
component score to the second regression model, 

[0174] Another preferred embodiment is a process for determining polymer 
properties in a fluidized-bed reactor system, the process comprising: (a) obtaining 
30 a locally weighted regression model for determining a polymer property selected 
from density, melt flow rate, molecular weight, molecular weight distribution, and 
functions thereof, the locally weighted regression model including principal 
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component loadings and principal component scores; (b) acquiring a Raman 
spectrum of a sample comprising polyolefin particles; (c) calculating a new 
principal component score from at least a portion of the Raman spectrum and the 
principal component loadings; and (d) calculating the polymer property by 
5 applying the new principal component score to the locally weighted regression 
model Even more preferred emodiments include one or more of the following: 
wherein the step of obtaining a regression model comprises: (i) obtaining a 
plurality of Raman spectra of samples comprising polyolefins; (ii) calculating 
principal component loadings and principal component scores from the spectra 

10 obtained in (i) using principal component analysis (PCA); and (iii) forming the 
regression model using the principal component scores calculated in (ii) such that 
the regression model correlates the polymer property to the principal component 
scores; wherein the step of acquiring a Raman spectrum comprises: (i) providing 
the sample of polyolefin particles; and (ii) irradiating the sample and collecting 

15 scattered radiation during a sampling interval using a sampling probe, wherein 
there is relative motion between the sample and the sampling probe during at least 
a portion of the sampling interval; wherein the process further comprises (i) 
obtaining a second regression model for determining a second polymer property, 
the second regression model including second principal component loadings and 
20 second principal component scores; (ii) calculating a new second principal 
component score from at least a portion of the Raman spectrum and the second 
principal component loadings; and (iii) calculating the second polymer property 
by applying the new second principal component score to the second regression 
model. 

25 [0175] Yet another preferred embodiment is a process for controlling polymer 
properties in a polymerization reactor system, the process comprising: (a) 
obtaining a regression model for determining a polymer property, the regression 
model including principal component loadings and principal component scores; 
(b) acquiring a Raman spectrum of a sample comprising polyolefin; (c) 

30 calculating a new principal component score from at least a portion of the Raman 
spectrum and the principal component loadings; (d) calculating the polymer 
property by applying the new principal component score to the regression model; 
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and (e) adjusting at least one polymerization parameter based on the calculated 
polymer property. Even more preferred emodiments include one or more of the 
following: wherein the step of obtaining a regression model comprises: (i) 
obtaining a plurality of Raman spectra of samples comprising polyolefins; (ii) 
calculating principal component loadings and principal component scores from 
the spectra obtained in (i) using principal component analysis (PCA); and (iii) 
forming the regression model using the principal component scores calculated in 
(ii) such that the regression model correlates the polymer property to the principal 
component scores; wherein the regression model is a locally weighted regression 
model; wherein the polymer property is selected from density, melt flow rate, 
molecular weight, molecular weight distribution, and functions thereof; wherein 
the sample comprises polyolefin particles; wherein the step of acquiring a Raman 
spectrum comprises: (i) providing the sample of polyolefin particles; and (ii) 
irradiating the sample and collecting scattered radiation during a sampling interval 
15 using a sampling probe, wherein there is relative motion between the sample and 
the sampling probe during at least a portion of the sampling interval; wherein the 
polymerization reactor is a fluidized-bed reactor; wherein the at least one 
polymerization parameter is. selected from the group consisting of monomer feed 
rate, comonomer feed rate, catalyst feed rate, hydrogen gas feed rate, and reaction 
temperature; the process further comprising: (i) obtaining a second regression 
model for determining a second polymer property, the second regression model 
including second principal component loadings and second principal component 
scores; (ii) calculating a new second principal component score from at least a 
portion of the Raman spectrum and the second principal component loadings; and 
25 (iii) calculating the second polymer property by applying the new second principal 
component score to the second regression model, and wherein the step of 
adjusting comprises adjusting at least one polymerization parameter based on the 
calculated polymer property, the calculated second polymer property, or both 
calculated polymer properties. 

10176] Yet still another preferred embodiment is a process for controlling polymer 
properties in a fluidized reactor system, the process comprising: (a) obtaining a 
locally weighted regression model for detemiining a polymer property selected 
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from density, melt flow rate, molecular weight, molecular weight distribution, and 
functions thereof, the locally weighted regression model including principal 
component loadings and principal component scores; (b) acquiring a Raman 
spectrum of a sample comprising polyolefin particles; (c) calculating a new 
principal component score from at least a portion of the Raman spectrum and the 
principal component loadings; (d) calculating the polymer property by applying 
the new principal component score to the locally weighted regression model; and 
(e) adjusting at least one polymerization parameter based on the calculated 
polymer property. Even more preferred emodiments include one or more of the 
following: wherein the step of obtaining a regression model comprises: (i) 
obtaining a plurality of Raman spectra of samples comprising polyolefins; (ii) 
calculating principal component loadings and principal component scores from 
the spectra obtained in (i) using principal component analysis (PCA); and (iii) 
forming the regression model using the principal component scores calculated in 
(ii) such that the regression model correlates the polymer property to the principal 
component scores; wherein the step of acquiring a Raman spectrum comprises: (i) 
providing the sample of polyolefin particles; and (ii) irradiating the sample and 
collecting scattered radiation during a sampling interval using a sampling probe, 
wherein there is relative motion between the sample and the sampling probe 
during at least a portion of the sampling interval; wherein the at least one 
polymerization parameter is selected from the group consisting of monomer feed 
rate, comonomer feed rate, catalyst feed rate, hydrogen gas feed rate, and reaction 
temperature; the process further comprising: (i) obtaining a second regression 
model for determining a second polymer property, the second regression model 
including second principal component loadings and second principal component 
scores; (ii) calculating a new second principal component score from at least a 
portion of the Raman spectrum and the second principal component loadings; and 
(iii) calculating the second polymer property by applying the new second principal 
component score to the second regression model, and wherein the step of 
adjusting comprises adjusting at least one polymerization parameter based on the 
calculated polymer property, the calculated second polymer property, or both 
calculated polymer properties. 
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[0177] An even more preferred embodiment of the invention includes any of the 
foregoing preferred embodiments, with or without the more preferred 
embodiments, wherein the Raman probe is inserted in situ into the polymerization 
reactor system, especially in a location where granular polymer is moving, for 
5 example inserted directly into the reactor body. Embodiments of this even more 
preferred embodiment include the following, either alone or in combination: 
wherein the polymerization reactor system is a gas phase polymerization reactor 
system; wherein the reactor body 22 is a fluidized bed reactor; wherein the Raman 
probe is purged with a stream of, for instance, N 2 or ethylene; wherein the 
10 aforementioned period of purging is cycled with a period of data collection; 
wherein the Raman probe is inserted in situ into at least one of the locations 
within the polymerization reactor system selected from the reactor body, the cycle 
gas piping, the product discharge system downstream of the reactor body, in the 
cyclone, in the purger/degasser, in the transfer line to finishing/pack-out, and in 

15 the feed bins to the extruder; and wherein the step of acquiring a Raman spectrum 
comprises: (ii) irradiating the sample of polymer, e.g., polyolefin, and collecting 
scattered radiation during a sampling interval using a Raman probe, and (ii) 
purging polymer from said Raman probe during a purging interval. 
[0178] A yet still more preferred embodiment of the foregoing even more 

20 preferred embodiment includes: (A) a gas phase polymerization reactor wherein 
gaseous monomer is introduced into a reactor body and polymer is discharged 
from the reactor, the improvement comprising a Raman probe inserted directly 
into said reactor body, whereby a Raman spectrum correlated to at least one 
polymer property is obtained; and (B) a gas phase polymerization process wherein 

25 gaseous monomer is introduced into a reactor body, and polymer is produced in 
said reactor body and polymer product is discharged from the reactor, the 
improvement comprising measuring at least one- property of the polymer produced 
in said reactor body by acquiring a Raman spectrum of said polymer within said 
reactor body. Yet even still more preferred embodiments of (B) include: wherein 

30 said Raman spectrum is acquired by inserting a Raman probe directly into said 
reactor body, and an optional probe purge, wherein said Raman probe is purged of 
polymer product by, for instance, a stream of nitrogen, ethylene (or monomer(s) 
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used in the polymerization reaction), hydrogen, and the like. The process also can 
include, among other variations that would be readily apparent to one of ordinary 
skill in the art with the present disclosure before them, (a) obtaining a regression 
model for determining a polymer property, the regression model including 
5 principal component loadings and principal component scores; (b) acquiring a 
Raman spectrum of a sample comprising polyolefin; (c) calculating a new 
principal component score from at least a portion of the Raman spectrum and the 
principal component loadings; and (d) calculating the polymer property by 
applying the new principal component score to the regression model; and further 

10 may comprise at least one polymerization parameter based on the polymer 
property, in moreover another very preferred embodiment wherein the at least one 
polymerization parameter is selected from at least one of the group consisting of 
monomer feed rate, comonomer (if present) feed rate, catalyst feed rate, hydrogen 
gas feed rate, reaction temperature. 

15 [0179] Possibly the most advantageous improvement provided by the present 
invention is illustrated by the following additional more preferred embodiment: 

(I) a gas phase polymerization reactor system wherein gaseous monomer is 
introduced into a reactor body and polymer is discharged from the reactor, the 
improvement comprising inserting a Raman probe in situ into said reactor system, 

20 whereby a Raman spectrum correlated to at least one property selected from the 
group consisting of a polymer property and a reactor operability property is 
obtained; including the embodiment wherein the Raman probe is inserted in situ 
into at least one of the locations within said polymerization reactor system 
selected from the group consisting of a polymerization reactor body, cycle gas 

25 piping, product discharge system downstream of the polymerization reactor body, 
a purger/degasser, a transfer line to finishing/pack-out, a feed bin to the extruder; 

(II) a gas phase polymerization process including a polymerization reactor system 
wherein gaseous monomer is introduced into a reactor body, polymer is produced 
in said reactor body, and polymer product is discharged from the reactor, the 

30 improvement comprising acquiring a Raman spectrum correlated with at least one 
property selected from the group consisting of a polymer property and a reactor 
operability property; and including the following embodiments, whose features 
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may be combined: wherein said Raman spectrum is acquired by a Raman probe 
inserted in situ into said polymerization reactor system, such as wherein the 
Raman probe is inserted in situ into at least one of the locations within said 
polymerization reactor system selected from the group consisting of a 
5 polymerization reactor body, cycle gas piping, product discharge system 
downstream of the polymerization reactor body, a purger/degasser, a transfer line 
to finishing/pack-out, a feed bin to the extruder; wherein the process further 
comprises purging polymer from said Raman probe, such as wherein said purging 
comprises purging with a stream of nitrogen gas; the process further comprising: 

10 (a) obtaining a regression model for determining a polymer property or a property 
correlated with reactor operability, the regression model including principal 
component loadings and principal component scores; (b) acquiring a Raman 
spectrum of a sample comprising polyolefin; (c) calculating a new principal 
component score from at least a portion of the Raman spectrum and the principal 

15 component loadings; and (d) calculating the polymer property or property 
correlated with reactor operability by applying the new principal component score 
to the regression model; any of the aforementioned further comprising adjusting at 
least one polymerization parameter based on the polymer property or property 
correlated with reactor operability, especially wherein the at least one 

20 polymerization parameter is selected from at least one of the group consisting of 
monomer feed rate, comonomer feed rate, catalyst feed rate, hydrogen gas feed 
rate, and reaction temperature. 

[0180] Preferred embodiments also include the apparatus including both the 
extractive sampling case as illustrated by Figure 2 and the in situ case described in 

25 detail above and includes: a gas phase polymerization reactor system, wherein 
gaseous monomer is introduced into a reactor body and polymer is discharged 
from the reactor, the improvement comprising providing a Raman probe in an 
extractive sampling system whereby a Raman spectrum correlated to at least one 
property selected from the group consisting of a polymer property and a reactor 

30 operability property is obtained, and more specifically wherein the extractive 
sampling system extracts polymer from a location selected from the group 
consisting of the cycle gas piping, the product discharge system downstream of 
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the exiting point of product, the transfer line between the product discharge 
system and the purger(s)/degasser(s), one or more of the purger(s)/degasser(s), the 
transfer line to finishing/pack-out, and the feed bins to the extruder/mixer; and 
also a gas phase polymerization reactor system wherein gaseous monomer is 
5 introduced into a reactor body and polymer is discharged from the reactor, the 
improvement comprising inserting a Raman probe in situ into said reactor system, 
whereby a Raman spectrum correlated to at least one property selected from the 
group consisting of a polymer property and a reactor operability property is 
obtained, and more specifically wherein the Raman probe is inserted in situ into at 

10 least one of the locations within said polymerization reactor system selected from 
the group consisting of the cycle gas piping, the product discharge system 
downstream of the exiting point of product, the transfer line between the product 
discharge system and the purger(s)/degasser(s), one or more of the 
purger(s)/degasser(s), the transfer line to finishing/pack-out, and the feed bins to 

15 the extruder/mixer. 

[0181] Finally, it should be noted that it may be particularly beneficial if the probe 
purge, if used, may be accomplished using a stream of nitrogen gas, monomer 
used in the polymerization reaction, or a combination of the two together or 
separately at different times and/or intervals. In addition it should be noted that in 

20 the extractive sampling technique described above, the sampling (e.g., as 
illustrated by Fig. 2) may be from one or more of the specific locations noted 
above for the in situ sampling case. 

[0182] Various tradenames used herein are indicated by a ™ symbol, indicating 
that the names may be protected by certain trademark rights. Some such names 
25 may also be registered trademarks in various jurisdictions. 

[0183] All patents, including the priority documents cited at the outset, and any 
other documents cited herein, such as ASTM or other test methods, are fully 
incorporated by reference to the extent such disclosure is not inconsistent with this 
invention and for all jurisdictions in which such incorporation is permitted. 
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CLAIMS 

What is claimed is: 

1. A process for determining polymer properties in a polymerization reactor 
system, the process comprising: 

(a) obtaining a regression model for determining a polymer property, 
the regression model including principal component loadings and 
principal component scores; 

(b) acquiring a Raman spectrum of a sample comprising polyolefin; 

(c) calculating a new principal component score from at least a portion 
of the Raman spectrum and the principal component loadings; and 

(d) calculating the polymer property by applying the new principal 
component score to the regression model; 

wherein the Raman spectrum acquired in step (b) is acquired from a 
Raman probe inserted in situ into the polymerization reactor system. 

2. The process of claim 1, wherein the step of obtaining a regression model 
comprises: 

(i) obtaining a plurality of Raman spectra of samples comprising 
polyolefins; 

(ii) calculating principal component loadings and principal component 
scores from the spectra obtained in (i) using principal component 
analysis (PC A); and 

(iii) forming the regression model using the principal component scores 
calculated in (ii) such that the regression model correlates the 
polymer property to the principal component scores. 

3. The process of claim 1, wherein the regression model is a locally weighted 
regression model. 

4. The process of claim 1, wherein the polymer property is selected from 
density, melt flow rate, molecular weight, molecular weight distribution, 
and functions thereof. 
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5. 
6. 

7. 
8. 



9. 



The process of claim 1, wherein the sample comprises polyolefin particles. 

The process of claim 1, wherein the step of acquiring a Raman spectrum 
comprises: (i) irradiating the sample of polyolefin and collecting scattered 
radiation during a sampling interval using a sampling probe, and (ii) 
purging polymer from said Raman .probe during a purging interval 

The process of claim 1, wherein the polymerization reactor is a fluidized- 
bed reactor. 

The process of claim 1, further comprising: 

(i) obtaining a second regression model for determining a second 
polymer property, the second regression model including second 
principal component loadings and second principal component 
scores; 

(ii) calculating a new second principal component score from at least a 
portion of the Raman spectrum and the second principal 
component loadings; and 

(iii) calculating the second polymer property by applying the new 
second principal component score to the second regression model. 

The process of claim 1, wherein the Raman probe is inserted in situ into 
said polymerization reactor system in a location where granular polymer is 



The process of claim 1, wherein the Raman probe is inserted in situ into at 
least one of the locations within said polymerization reactor system 
selected from the group consisting of the cycle gas piping, the product 
discharge system downstream of the exiting point of product, the transfer 
line between the product discharge system and the purger(s)/degasser(s), 



moving. 
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one or more of the purger(s)/degasser(s), the transfer line to 
finishing/pack-out, and the feed bins to the extruder/mixer. 

11. The process of claim 1, wherein the Raman probe is inserted in situ into 
the reactor body. 

12. The process of claim 1, further comprising purging polymer from said 
Raman probe. 

13. A process for determining polymer properties in a fluidized-bed reactor 
system, the process comprising: 

(a) obtaining a locally weighted regression model for determining a 
polymer property selected from density, melt flow rate, molecular 
weight, molecular weight distribution, and functions thereof, the 
locally weighted regression model including principal component 
loadings and principal component scores; 

(b) acquiring a Raman spectrum of a sample comprising polyolefin 
particles; 

(c) calculating a new principal component score from at least a portion 
of the Raman spectrum and the principal component loadings; and 

(d) calculating the polymer property by applying the new principal 
component score to the locally weighted regression model 

wherein the Raman spectrum acquired in step (b) is acquired from a 
Raman probe inserted in situ into the polymerization reactor system. 

14. The process of claim 13, wherein the step of obtaining a regression model 
comprises: 

(i) obtaining a plurality of Raman spectra of samples comprising 
polyolefins; 

(ii) calculating principal component loadings and principal component 
scores from the spectra obtained in (i) using principal component 
analysis (PCA); and 
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(iii) fonning the regression model using the principal component scores 
calculated in (ii) such that the regression model correlates the 
polymer property to the principal component scores. 
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15. The process of claim 13, wherein the step of acquiring a Raman spectrum 
comprises: 

(i) providing the sample of polyolefin particles; and 

(ii) irradiating the sample and collecting scattered radiation during a 
sampling interval using a sampling probe, 

wherein there is relative motion between the sample and the sampling 
probe during at least a portion of the sampling interval. 



1 6. The process of claim 1 3, further comprising: 

(i) obtaining a second regression model for determining a second 
polymer property, the second regression model including second 
principal component loadings and second principal component 
scores; 

(ii) calculating a new second principal component score from at least a 
portion of the Raman spectrum and the second principal 
component loadings; and 

(iii) calculating the second polymer property by applying the new 
second principal component score to the second regression model. 

17. The process of claim 13, wherein the Raman probe is inserted in situ into 
said polymerization reactor system in a location where granular polymer is 
moving. 

18. The process of claim 13, wherein the Raman probe is inserted in situ into, 
at least one of the locations within said polymerization reactor system 
selected from the group consisting of the cycle gas piping, the product 
discharge system downstream of the exiting point of product, the transfer 
line between the product discharge system and the purger(s)/degasser(s), 
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one or more of the purger(s)/degasser(s), the transfer line to 
finishing/pack-out, and the feed bins to the extruder/mixer. 

19. The process of claim 13, wherein the Raman 1 probe is inserted in situ into 
the reactor body. 

20. The process of claim 13, further comprising a step of purging polymer 
from said Raman probe. 

21. A process for controlling polymer properties in a polymerization reactor 
system, the process comprising: 

(a) obtaining a regression model for determining a polymer property, 
the regression model including principal component loadings and 
principal component scores; 

(b) acquiring a Raman spectrum of a sample comprising polyolefin; 

(c) calculating a new principal component score from at least a portion 
of the Raman spectrum and the principal component loadings; 

(d) calculating the polymer property by applying the new principal 
component score to the regression model; and 

(e) adjusting at least one polymerization parameter based on the 
calculated polymer property; 

wherein the Raman spectrum acquired in step (b) is acquired from a 
Raman probe inserted in situ into the polymerization reactor system. 

22. The process of claim 21, wherein the step of obtaining a regression model 
comprises: 

(i) obtaining a plurality of Raman spectra of samples comprising 
polyolefins; 

(ii) calculating principal component loadings and principal component 
scores from the spectra obtained in (i) using principal component 
analysis (PCA); and 
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(iii) forming the regression model using the principal component scores 
calculated in (ii) such that the regression model correlates the 
polymer property to the principal component scores. 

23. The process of claim 21, wherein the regression model is a locally 
weighted regression model. 

24. The process of claim 21, wherein the polymer property is selected from 
density, melt flow rate, molecular weight, molecular weight distribution, 
and Junctions thereof. 

25. The process of claim 21, wherein the sample comprises polyolefin 
particles. 

26. The process of claim 21, wherein the step of acquiring a Raman spectrum 
comprises: (i) irradiating the sample of polyolefin and collecting scatteredj 
radiation during a sampling interval using a sampling probe, and (ii) 
purging polymer from said Raman probe during a purging interval. 

27. The process of claim 21, wherein the polymerization reactor is a fluidized- 
bed reactor. 

28. The process of claim 21, wherein the at least one polymerization parameter 
is selected from the group consisting of monomer feed rate, comonomer 
feed rate, catalyst feed rate, hydrogen gas feed rate, and reaction 
temperature. 

29. The process of claim 2 1 , further comprising: 

(i) obtaining a second regression model for determining a second 
polymer property, the second regression model including second 
principal component loadings and second principal component 
scores; 
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(ii) calculating a new second principal component score from at least a 
portion of the Raman spectrum and the second principal 



(iii) calculating the second polymer property by applying the new 

second principal component score to the second regression model, 
and wherein the step of adjusting comprises adjusting at least one 
polymerization parameter based on the calculated polymer property, the 
calculated second polymer property, or both calculated polymer properties. 

A process for controlling polymer properties in a fluidized reactor system, 
the process comprising: 

(a) obtaining a locally weighted regression model for determining a 
polymer property selected from density, melt flow rate, molecular 
weight, molecular weight distribution, and functions thereof, the 
locally weighted regression model including principal component 
loadings and principal component scores; 

(b) acquiring a Raman spectrum of a sample comprising polyolefin 
particles; 

(c) calculating a new principal component score from at least a portion 
of the Raman spectrum and the principal component loadings; 

(d) calculating the polymer property by applying the new principal 
component score to the locally weighted regression model; and 

(e) adjusting at least one polymerization parameter based on the 
calculated polymer property; 

wherein the Raman spectrum acquired in step (b) is acquired from a 
Raman probe inserted in situ into the polymerization reactor system. 

The process of claim 30, wherein the step of obtaining a regression model 
comprises: 

(i) obtaining a plurality of Raman spectra of samples comprising 
polyolefins; 



component loadings; and 
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(ii) calculating principal component loadings and principal component 
scores from the spectra obtained in (i) using principal component 
analysis (PCA); and 

(iii) forming the regression model using the principal component scores 
calculated in (ii) such that the regression model correlates the 
polymer property to the principal component scores. 

32. The process of claim 30, wherein the step of acquiring a Raman spectrum 
comprises: 

(i) providing the sample of polyolefin particles; and 

(ii) irradiating the sample and collecting scattered radiation during a 
sampling interval using a sampling probe, 

wherein there is relative motion between the sample and the sampling 
probe during at least a portion of the sampling interval. 

33. The process of claim 30, wherein the at least one polymerization parameter 
is selected from the group consisting of monomer feed rate, comonomer 
feed rate, catalyst feed rate, hydrogen gas feed rate, and reaction 
temperature. 

34 The process of claim 30, further comprising: 

(i) obtaining a second regression model for determining a second 
polymer property, the second regression model including second 
principal component loadings and second principal component 
scores; 

(ii) calculating a new second principal component score from at least a 
portion of the Raman spectrum and the second principal 
component loadings; and 

(iii) calculating the second polymer property by applying the new 
second principal component score to the second regression model, 
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and wherein the step of adjusting comprises adjusting at least one 
polymerization parameter based on the calculated polymer property, the 
calculated second polymer property, or both calculated polymer properties. 

35. In a gas phase polymerization reactor system wherein gaseous monomer is 
introduced into a reactor body and polymer is discharged from the reactor, 
the improvement comprising inserting a Raman probe in situ into said 
reactor system, whereby a Raman spectrum correlated to at least one 
property selected from the group consisting of a polymer property and a 
reactor operability property is obtained. 

i 

36. The gas phase polymerization reactor system according to Claim 35, 
wherein the Raman probe is inserted in situ into at least one of the 
locations within said polymerization reactor system selected from the 
group consisting of the cycle gas piping, the product discharge system 
downstream of the exiting point of product, the transfer line between the 
product discharge system and the purger(s)/degasser(s), one or more of the 
purger(s)/degasser(s), the transfer line to finishing/pack-out, and the feed 
bins to the extruder/mixer. 

37. In a gas phase polymerization process including a polymerization reactor 
system wherein gaseous monomer is introduced into a reactor body, 
polymer is produced in said reactor body, and polymer product is 
discharged from the reactor, the improvement comprising acquiring a 
Raman spectrum correlated with at least one property selected from the 
group consisting of a polymer property and a reactor operability property. 

38. The process according to claim 37, wherein said Raman spectrum is 
acquired by a Raman probe inserted in situ into said polymerization 
reactor system. 
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39. The process according to claim 38, wherein the Raman probe is inserted in 
situ into at least one of the locations within said polymerization reactor 
system selected from the group consisting of the cycle gas piping, the 
product discharge system downstream of the exiting point of product, the 
transfer line between the product discharge system and the 
purger(s)/degasser(s), one or more of the purger(s)/degasser(s), the transfer 
line to finishing/pack-out, and the feed bins to the extruder/mixer.. 

40. The process according to claim 38, further comprising purging polymer 
from said Raman probe. 

41. The process according to claim 38, wherein said purging comprises 
purging with a stream of nitrogen gas or monomer. 

42. The process according to claim 38, further comprising: 

(a) obtaining a regression model for determining a polymer property or 
a property correlated with reactor operability, the regression model 
including principal component loadings and principal component 
scores; 

(b) acquiring a Raman spectrum of a sample comprising polyolefin; 

(c) calculating a new principal component score from at least a portion 
of the Raman spectrum and the principal component loadings; and 

(d) calculating the polymer property or property correlated with reactor 
operability by applying the new principal component score to the 
regression model. 

43. The process according to claim 42, further comprising adjusting at least 
one polymerization parameter based on the polymer property or property 
correlated with reactor operability. 



44. 



The process according to claim 43, wherein the at least one polymerization 
parameter is selected from at least one of the group consisting of monomer 
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feed rate, comonomer feed rate, catalyst feed rate, hydrogen gas feed rate, 
and reaction temperature. 

45. The process according to claim 37, wherein said Raman spectrum is 
acquired by extractive sampling from said polymerization reactor system. 



46. 



The process according to claim 45, wherein said extractive sampling is 
from at least one of the locations within said polymerization reactor 
system selected from the group consisting of the cycle gas piping, the 
product discharge system downstream of the exiting point of product, the 
transfer line between the product discharge system and the 
purger(s)/degasser(s), one or more of the purger(s)/degasser(s), the transfer 
line to finishing/pack-out, and the feed bins to the extruder/mixer. 



is 



47. In a gas phase polymerization reactor system wherein gaseous monomer 
introduced into a reactor body and polymer is discharged from the reactor, 
the improvement comprising providing a Raman probe in an extractive 
sampling system whereby a Raman spectrum correlated to at least one 
property selected from the group consisting of a polymer property and a 
reactor operability property is obtained. 



The gas phase polymerization reactor system according to Claim 47, 
wherein the extractive sampling system extracts polymer from a location 
selected from the group consisting of the cycle gas piping, the product 
discharge system downstream of the exiting point of product, the transfer 
line between the product discharge system and the purger(s)/degasser(s), 
one or more of the purger(s)/degasser(s), the transfer line to 
finishing/pack-out, and the feed bins to the extruder/mixer. 
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